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ABSTRACT 

The  objective  of  this  study  is  to  develop  classification 
equations  to  forecast  the  daily  probability  of  occurrence  of 
marine  fog  at  selected  locations  on  the  west  coast  of  the 
United  States,  using  parameters  easily  obtainable  from  upper- 
air  soundings  and  surface  observations.   In  order  to  achieve 
this  objective  a  computerized  stepv;ise  linear  discriminant 
analysis  program  is  extensively  employed.   Data  input  con- 
sists of  surface  and  radiosonde  observations  for  the  five- 
year  period  1  July  1968  to  30  June  1973  at  three  U.  S.  west 
coast  stations,  namely  San  Diego  and  Oakland,  California  and 
Quillayute,  Washington. 

Tables  showing  the  number  of  fog  and  no-fog  cases,  the 
classification  functions,  and  the  percentages  of  correct  fog 
and  no-fog  discrimination  are  presented  for  each  station. 
The  most  capable  fog/no-fog  discrimination  parameters  are 
discussed  for  each  set  of  classification  equations.   Test 
results  for  the  San  Diego  equations  using  a  three-year  inde- 
pendent data  set  are  also  shown. 
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I.    INTRODUCTION 

Marine  fog  is  an  impediment  to  both  commercial  and  mili- 
tary shipping  and  air  operations.   Such  fog  has  often  caused 
costly  delays  to  commercial  ship  owners  and  domestic  air 
carriers  while  endangering  lives  and  property  as  a  result  of 
low  visibilities.   Naval  operations  such  as  aircraft  transfer 
and  landing,  antisubmarine  warfare  and  amphibious  maneuvers, 
navigational  reconnaissance  and  search  and  rescue  missions 
may  be  hampered  by  marine  fog.   The  importance  and  impact  of 
marine  fog  on  naval  operations  during  World  War  II  and  the 
recent  costs  (1969-1974)  in  lives  and  property  damage  of 
ship  and  aircraft  accidents  suffered  by  the  United  States 
Navy  during  fog  situations  has  been  described  by  Wheeler 
(1974)  . 

Recognized  as  an  important  meteorological  problem,  marine 
fog  prediction  and  dissipation  have  been  the  subjects  of  on- 
going research  by  several  Department  of  Defense  (DOD)  acti- 
vities.  In  particular  the  Naval  Air  Systems  Command  (NAVAIR) 
and  the  Office  of  Naval  Research  (ONR)  have  jointly  sponsored 
marine  fog  research  in  an  attempt  to  model  and  simulate  the 
fog  for  purposes  of  prediction.   A  Naval  Postgradute  School 
(NPS)  group,  under  the  direction  of  Dr.  Leipper,  Chairman  of 
the  Department  of  Oceanography,  and  Dr.  Renard  of  the  Depart- 
ment of  Meteorology,  is  participating  in  this  coordinated 
effort  with  the  following  specific  objectives: 
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"1.   To  observe  and  describe  the  formation  and  dissipa- 
tion of  fog  at  sea  in  cooperation  with  Cornell 
Aeronautical  Laboratory,  Inc. 
2.   To  relate  the  synoptic  variables  in  meteorology 

and  oceanography  which  are  important  in  marine  fog 
formation  and  dissipation  to  microphy sical  and 
micrometeorological  characteristics . 
3-   To  improve  marine  fog  forecasting  and  to  develop 
marine  fog  prediction  methods." 
Several  recent  NPS  studies  have  addressed  the  problem  of 
marine  fog  climatology  of  the  North  Pacific  Ocean  and  incor- 
porated a  method  of  synthesizing  ships'  surface  synoptic  re- 
ports for  the  purpose  of  deriving  frequencies  of  marine  fog 
occurrence  (Englebr e t son ,  1974;  Daughenbaugh ,  1975;  Renard, 
Englebretson  and  Daughenbaugh,  1975;  Willms,  1975)  .   Wallace 
(1975)  and  Hale  (1975)  have  explored  the  feasibility  of  using 
meteorological  satellite  data  as  a  means  of  discerning  the 
presence  of  marine  fog.   In  a  study  by  McClure  (1974) ,  a  po- 
tential link  between  temperature  inversions  and  the  occur- 
rence of  marine  advection  fog  was  investigated,  in  order  that 
temperature  inversion  parameters  could  possibly  be  used  as 
forecasting  aids  at  coastal  stations. 

Although  the  Environmental  Prediction  Research  Facility 
(EPRF),  Monterey,  California,  is  engaged  in  numerical 


Leipper,  D.  F.,  and  others,  1973:  Observation,  Analysis 
and  Prediction  of  Marine  Fog.   Naval  Postgraduate  School, 
Department  of  Oceanography  Annual  Report  for  period  1  July 
1972  to  30  June  1973,  p.  10. 
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modeling  studies  on  the  physics  of  fog  formation,  mainte- 
nance and  dissipation  (Feit,  1972;  Barker,  1973),  and 
researchers  elsewhere,  such  as  Mack  et  al  (1973)  ,  have 
developed  numerical  models  of  advection  fog,  few  opera- 
tional forecasting  schemes  for  marine  fog  on  the  west 
coast  of  the  United  States  have  been  developed  to  date. 
Fleet  Numerical  Weather  Central  (FNWC),  Monterey,  Cali- 
fornia, has  an  operational  advection  fog  program  ("FTER") 
which  forecasts  grid  point  fog  probabilities  to  48  hours 
based  on  a  system  of  multi-parameter  tests  (U.  S.  Naval 
Weather  Service  Command,  1975) .   Given  the  correct  humid- 
ity values,  this  program  gives  acceptable  results  in  cases 
of  broad-scale  fog.   However,  its  primary  weakness  stems 
from  the  large  63x63  grid  --  local  effects  just  cannot  be 
reconciled . 

"Forecasting  fog  is  one  of  meteorology's  hardest  prob- 
lems.  This  is  because  there  are  so  many  factors  involved-- 

moisture,  temperature,  wind,  geography,  stability,  cloud 

2 
cover,  and  several  more."    It  is  hoped  that  this  study 

may  provide  insight  and  methodology  which  will  aid  fore- 
casters in  predicting  the  occurrence  of  marine  fog  for 
selected  locations  on  the  west  coast  of  the  United  States. 


U.  S.  Naval  Weather  Service  Command,  1975:   U.  S. 
Naval  Weather  Service  Numerical  Environmental  Products 

Manual,  NAVAIR  ,  Department  of  the  Navy, 

Washington,  D.  C.,  p.  4.26-2. 
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II.   BACKGROUND 

A,   A  REVIEW  OF  SELECTED  PREVIOUS  FOG  RESEARCH 

A  considerable  amount  of  research  into  the  problem  of 
marine  fog  on  the  North  American  west  coast  has  been  con- 
ducted in  the  southern  California  area.   In  1948  Leipper 
(1948)  devised  a  model  to  illustrate  the  manner  in  which 
fog  situations  develop  in  the  San  Diego  area.   This  model 
was  described  in  four  stages  and  is  still  used  to  some  de- 
gree by  military  forecasters  in  the  southern  California 
area  . 

The  first  stage  or  initial  conditions  in  Leipper ' s  model 
requires  the  presence  of  air  over  the  sea  which  has  a  tem- 
perature higher  than  the  sea-surface  temperature  and  which 
is  relatively  dry  aloft.   This  condition  guarantees  the  for- 
mation of  an  inversion  which  will  restrict  the  vertical 
movement  of  moisture  and  thus  cause  a  thin  lower  layer  of 
air  to  approach  saturation.   Such  initial  conditions  are 
typically  brought  about  when  a  lobe  of  the  North  Pacific  sub- 
tropical anticyclone  pushes  inland  over  northern  California 
causing  a  general  easterly  flow  over  southern  California  so 
that  the  air  arriving  at  San  Diego  is  warm  and  dry,  having 
descended  the  4000-foot  slope  of  the  coastal  mountains. 

In  the  second  stage  of  fog  development,  the  easterly 
flow  weakens,  the  flow  of  warm  subsiding  air  toward  the  sea 
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decreases,  and  the  air  offshore  from  San  Diego  remains  in  a 
relatively  stagnant  condition.   Due  to  evaporation  from  the 
sea  the  moisture  content  of  the  lower  layers  increases  con- 
siderably, and  the  heat  content  decreases  due  to  the  conduc- 
tion of  heat  downward.   During  this  process  a  surface 
inversion  is  formed  and  the  lowest  air  layers  become  nearly 
saturated  at  a  temperature  close  to  that  of  the  underlying 
sea  sur f ac e . 

Fog  forms  in  the  third  stage  as  the  normal  northwesterly 
airflov/  and  sea  breeze  regime  returns.   The  fog  is  created 
in  the  thin  surface  layer  at  sea  as  the  v;arm,  nearly  satu- 
rated air  moves  over  the  colder  coastal  waters.   Thus,  fog 
is  produced  by  cooling  from  below  (Leipper,  1968) .   Once  the 
white  cloud  is  formed,  incoming  radiation  reflected  from  the 
top  and  heat  radiated  from  the  cloud  itself  results  in  cool- 
ing of  the  thin  fog  layer  to  a  temperature  lower  than  the 
ocean  surface  in  the  vicinity.   Then  as  the  sea  breeze 
strengthens  during  the  daytime,  the  fog  is  advected  shore- 
ward . 

The  fourth  and  final  stage  is  characterized  by  the  deepen- 
ing of  the  nearly  adiabatic  lower  layer  to  approximately 
400  feet  or  more.   Since  evaporation  from  the  sea  continues 
and  the  presence  of  the  strong  inversion   still  restricts 
the  vertical  movement  of  moisture,  the  dew  point  remains  high 
in  the  marine  layer.   After  the  daytime  sea  breeze  carries 
the  fog  over  the  shore,  radiation  and  conduction  of  heat  from 
the  land  mass  dissipate  the  fog  somewhat,  but  the  nearly 
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adiabatic  layer,  being  deeper  than  the  previous  stage,  is 
not  destroyed  as  rapidly.   Rather,  it  retains  its  high 
moisture  content  in  vapor  form  and  moves  inland.   As  the 
marine  layer  beneath  the  inversion  continues  to  deepen,  a 
point  is  reached  where  the  mixing  and  cooling  processes  do 
not  result  in  condensation  throughout  the  layer.   Then  the 
upper  portion  fills  with  water  particles  first  and  fog 
does  not  form  on  the  ground.   In  this  way  the  fog  sequence 
is  ended  and  the  stratus  regime  begins^   The  complete  fog 
sequence  Leipper  described  usually  extends  over  a  period  of 
about  five  days. 

In  a  later  work  Leipper  (1968)  compared  the  meteorologi- 
cal conditions  associated  with  the  observance  of  a  sharp 
smog  bank,  near  Riverside,  California  with  those  previously 
shown  to  be  related  to  winter  fogs  in  southern  California. 
He  believed  that  the  characteristics  of  certain  stages  of 
these  fog  situations  in  winter  were  quite  similar  to  those 
observed  by  Edinger  (1963)  in  July  1961  and  to  those  re- 
ported by  Stephens  (1965).   The  distinguishing  characteris- 
tic was  the  presence  of  unusually  warm,  dry  air  over  the 
inversion  and  a  very  cold  marine  layer  beneath  it.   In  such 
situations  the  air  above  the  inversion  is  often  nearly  10  C 
warmer  than  the  underlying  sea.   The  combination  of  very 
warm  air  aloft  and  an  unusually  cold  marine  layer  results 
in  the  strongly  stable  situation  described  by  both  Edinger 
and  Stephens  as  having  high  smog  potential  when  it  moves  in- 
land.  Such  situations  in  the  San  Diego  area  were  associated 
with  fogs  having  sharp  boundaries. 
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Three  nondiurnal  indices  presented  by  Leipper  in  1948, 
and  found  useful  in  the  prediction  of  west  coast  fog,  also 
may  be  important  in  predicting  situations  favorable  to  the 
shallow  sharp-banked  smogs  which  are  often  observed.   The 
indices  are  composed  from  data  obtained  in  the  morning  radio- 
sondes used  together  with  the  coastal  water  temperature  and 
the  surface  dew  point  in  the  marine  layer.   Leipper  (1948) 
described  the  indices  as  follows: 

1.  Height  of  the  Inversion  Base:   The  height  above 
which  the  air  temperature  increases  with  height  at 
the  most  rapid  rate  on  the  morning  raob  (radiosonde 
observation) . 

2.  Temperature  Index:   The  quantity  (Ta-Tw);  if  an  in- 
version exists  with  base  below  3,000  feet,  Ta  is  the 
highest  air  temperature  above  the  inversion  base  on 
the  morning  raob,  otherwise  Ta  is  the  surface  air 
temperature,  and  Tw  is  the  coastal  sea-surface  tem- 
perature . 

3.  Moisture  Index:   The  difference  between  the  afternoon 
dew-point  and  the  coastal  water  temperatures. 

The  conditions  found  most  favorable  for  fog  occurrence, 
if  applicable  to  smog,  require  the  base  of  the  inversion  to 
be  below  1,300  feet,  the  temperature  index  to  be  positive, 
and  the  moisture  index  to  be  greater  than  -5  C, 

Schroeder  et  al  (1967)  describe  three  significant  aspects 
of  marine  air  invasion  along  the  west  coast  of  the  United 
States,  two  of  which  are  related  to  fog  development.   The 
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most  noticeable  invasion  is  that  of  the  diurnal  sea  breeze. 
This  well  known  phenomenon  is  produced  by  differential  heat- 
ing of  the  land  and  water  masses  and  results  in  the  sea 
breeze  circulation  cell  which  brings  cool,  moist  marine  air 
landward.   One  of  the  most  distinctive  characteristics  of 
the  sea  breeze  is  the  front-like  appearance  at  its  leading 
edge  . 

The  inland  movement  of  the  sea  breeze  and  its  front  is 
strongly  controlled  by  terrain  features.   In  addition,  the 
large-scale  weather  patterns  affect  the  sea  breeze.   When  the 
synoptic-scale  gradient  flow  is  directed  offshore,  the  sea- 
breeze  front  is  intense  but  does  not  penetrate  far  inland, 
and  the  strength  of  the  flow  is  weak.   However,  when  the 
gradient  flow  is  onshore,  the  front  is  weak  and  the  sea- 
breeze  flow  is  strong,  resulting  in  considerable  inland 
penetration . 

Perhaps  the  most  significant  marine  air  invasion  is  that 
of  the  United  States  west  coast  monsoon.   This  phenomenon  be- 
gins in  the  late  spring  and  continues  until  fall,  bringing  a 
slow,  steady  transport  of  marine  air  inland  from  the  North 
Pacific  Ocean  high  pressure  cell.   This  subtropical  high  is 
composed  of  a  shallow  marine  layer  capped  by  a  subsidence 
inversion.   The  marine  inversion  is  the  most  pronounced  and 
shallow  during  the  monsoon  season.   Irregular  fluctuations 
in  the  marine  inversion  are  generally  related  to  the  synoptic 
weather  patterns,  while  the  diurnal  variations  of  the  inver- 
sion depth  are  partly  related  to  the  interaction  of  the  sea 
breeze  and  the  monsoon.   The  monsoon  is  generally  confined  to 
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a  layer  less  than  2  km  deep  and  its  inland  path  is  primarily 
through  low  gaps  in  the  coastal  mountains. 

In  summary,  the  monsoon  is  a  feature  of  the  general 
atmospheric  circulation  which  undergoes  modification  at  the 
coastline,  interacts  with  both  the  sea  breeze  and  the  moving 
synoptic-scale  systems,  and  is  confined  to  a  very  shallow 
layer  o 

One  of  the  most  comprehensive  studies  of  marine-fog  fore- 
casting on  the  west  coast  was  done  at  Scripps  Institution 
of  Oceanography  during  a  contract  with  ONR,  and  was  super- 
vised by  D.  F.  Leipper  (Leipper  and  others,  1948) .   The 
purpose  of  this  "Fog  Project"  was  to  develop  and  to  test 
principles  which  might  serve  as  a  basis  for  fog  forecasting 
in  coastal  areas  and  possibly  elsewhere. 

The  "Final  Report  of  the  Fog  Project"  outlines  some  of 
the  theory  of  condensation  and  evaporation  which  have  parti- 
cular importance  in  the  development  and  application  of  the 
fog  forecasting  technique.   It  is  generally  believed  that 
the  local  fog  forecasting  problem  is  primarily  one  of  pre- 
dicting the  moisture  content  of  the  air  and  the  amount  of 
local  cooling  which  will  occur.   The  sea-surface  temperature 
is  an  important  variable  in  determining  the  amount  of  modi- 
fication (measured  by  the  change  in  content  of  heat  and 
moisture)  which  will  be  brought  about  in  a  given  overlying 
air  mass.   This  temperature  is  also  a  convenient  quantity 
with  which  other  more  erratic  quantities  such  as  air  tempera- 
ture and  dew  point  may  be  compared,  as  was  done  in  Leipper's 
(1948)  non-diurnal  fog  indices.   In  low  coastal  regions  the 
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dew-point  temperature  is  a  reliable  index  of  moisture  con- 
tent, and  is  the  most  convenient  quantity  for  use  in  fore- 
casting changes  in  moisture  which  will  affect  fog  formation. 

The  forecasting  procedure  devised  by  the  researchers  of 
the  "Fog  Project"  consists  of  three  parts.   First  the  fore- 
caster must  prepare  prognostic  charts  in  order  to  determine 
from  them  certain  fundamental  weather  features  (local  fea- 
tures such  as  wind  speed,  cloud  cover,  range  of  ceiling  and 
surface  pressure;  also,  synoptic-scale  features  such  as 
changes  in  surface  air  flow  and  subsidence) .   The  second  step 
involves  applying  statistical  aids  (e.g.,  c lima tological  data, 
temperature  and  time  forecasts,  and  visibility  data)  to  deter- 
mine the  local  changes  in  heat  and  moisture  content  of  the 
air  associated  with  the  fundamental  features  which  are  fore- 
cast.  Finally,  the  forecaster  must  determine  the  visibility 
from  the  moisture  content,  changes  in  heat  content  and  other 
available  data  by  applying  empirical  rules. 

Information  obtained  from  radiosondes  such  as  the  height 
of  the  inversion  base,  the  strength  and  thickness  of  the  in- 
version, etc.,  may  be  significant  in  fog  forecasting.   In  the 
"Final  Project"  J.  B.  Wickham  discusses  why  the  inversion  in- 
dex (height  of  the  inversion  base)  is  important.   Generally 
speaking,  when  a  strong  temperature  inversion  exists  over  a 
nearly  adiabatic  layer,  the  moisture  flux  through  the  inver- 
sion is  small  as  compared  with  the  flux  through  the  surface 
layer.   Consequently,  moisture  is  trapped  below  the  inversion 
and  over  a  period  of  time,  near  saturation  of  the  entire 
adiabatic  layer  can  occur.   Not  only  does  near  saturation 
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occur  most  rapidly  with  low  inversions,  but  the  amount  of 
cooling  necessary  to  produce  condensation  throughout  the 
layer  is  less  for  thin  layers  than  for  thick  ones.   Thus, 
the  time  required  to  approach  saturation  and  the  amount  of 
cooling  needed  to  form  fog  are  both  functions  of  the  height 
of  the  inversion  base.   The  inversion  index  is  of  further 
importance  since  it  effects  the  diurnal  temperature  range 
and  indicates  the  location  of  the  dry  marine-air  boundary. 
A  more  recent  investigation  into  the  relationships  be- 
tv/een  temperature  inversions  and  the  occurrence  of  marine 
advection  fog  was  undertaken  by  McClure  (1974) .   He  studied 
three  primary  inversion  parameters:   the  height  of  the  in- 
version base,  thickness  of  the  inversion  layer,  and  the  tem- 
perature gradient  within  the  inversion  layer,  using  ten 
months  of  surface  and  upper-air  data  (July  1973  to  April 
1974) .   The  upper-air  sounding  stations  used  were  San  Diego 

(Montgomery  Field),  California;  Oakland  (International  Air- 
port), California;  and  Quillayute,  Washington.   The  surface 
observation  stations  were  San  Diego  (Lindberg  Field),  Cal- 
ifornia; Oakland  (International  Airport),  California;  and 
Seattle-Tacoma  (International  Airport),  Washington. 

Because  of  the  considerable  geographic  separation 

(nearly  100  miles)  between  the  upper-air  soundings  taken  at 
Quillayute  and  the  surface  fog  observed  at  the  Seattle-Tacoma 
Airport,  no  definitive  relationships  v/ere  achieved.   However, 
at  San  Diego  a  relationship  between  the  change  in  inversion 
parameters  and  the  occurrence  of  fog  was  found  to  exist  in 
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both  the  summer  and  winter  seasons.   At  Oakland,  consistency 
in  the  change  of  one  parameter  was  shown  for  most  of  the 
data  period. 

A  brief  summary  of  some  of  McClure's  relationships  between 
temperature  inversions  and  marine  advection  fog  follows. 

1.  Summer  fog  cases  are  generally  preceded  by  an  inver- 
sion in  the  atmosphere  below  850  mb .   For  soundings 
taken  during  fog,  the  lack  of  an  inversion  below 
850  mb  is  the  exception  rather  than  the  rule, 
especially  during  the  summer  and  early  fall. 

2.  Most  fall  and  winter  fog  cases  (80%-90%)  are  pre- 
ceded by  an  inversion. 

3.  During  the  summer  season  in  San  Diego  the  thickness 
of  the  0000  Greenwich  Mean  Time  (GMT)  inversion 
layer  usually  decreases  within  24  hours  prior  to  the 
occurrence  of  fog  and  the  inversion  base  lowers. 
Cooling  of  the  surface  layer  precedes  the  fog.  Where 
fog  persists  there  is  a  general  warming  trend  in  the 
lower  atmosphere  associated  with  a  subsidence  inver- 
sion in  the  0000  GMT  sounding. 

4.  During  the  winter  season  in  San  Diego  the  inversion 
base  rises  in  the  last  24-hour  change  of  the  0000 
GMT  sounding  before  a  fog  sequence.   In  the  1200  GMT 
sounding  the  inversion  is  based  at  the  surface  for 
at  least  four  soundings  prior  to  the  fog.   Within 

24  hours  of  the  commencement  of  fog,  the  strength  of 
the  inversion  layer  lessens,  the  depth  increases, 
while  the  temperature  gradient  decreases. 
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5.  In  Oakland's  summer  season  there  is  a  decrease  in 
the  inversion  thickness  and  an  increase  in  the  tem- 
perature gradient  in  the  0000  GMT  sounding  within 
24  hours  prior  to  the  fog  occurrence. 

6.  In  Oakland's  fall  season  there  is  almost  always  a 
surface  inversion  in  the  last  two  1200  GMT  soundings 
prior  to  the  onset  of  fog.   This  is  also  true  for 
the  winter  months  of  January  and  February. 

7.  During  the  winter  season  in  Oakland  inversions  do 
not  always  accompany  or  precede  fog.   If  an  inversion 
exists,  it  is  usually  weakening  prior  to  fog  occur- 
rence, while  the  inversion  thickness  and  temperature 
gradient  are  decreasing. 

McClure's  investigations  and  recommendations  served  as 
the  platform  from  which  this  study  was  launched.   A  much 
broader  data  base  and  computer  analysis  was  believed  to  be 
essential  if  more  concrete  relationships  were  to  be  made. 
Additional  parameters  were  needed,  especially  within  the 
surface  observations.   Finally,  it  was  considered  to  be 
highly  desirable  to  have  the  upper-air  sounding  station  and 
the  surface  observation  station  at  the  same  location.   These 
changes  and  others,  plus  a  different  analysis  approach,  were 
incorporated  in  this  research. 
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Ill .   OBJECTIVE,   DATA  AND  DEFINITIONS 

A.  OBJECTIVE  AND  APPROACH 

The  objective  of  this  study  was  to  develop  classifica- 
tion equations  to  forecast  the  daily  probability  of  occur- 
rence of  marine  fog  at  selected  locations  on  the  west  coast 
of  the  United  States,  using  parameters  easily  obtainable 
from  upper-air  soundings  and  surface  observations.   In  order 
to  achieve  this  objective,  a  stepwise  linear  discriminant 
analysis  program,  BMD07M  (Dixon,  1973),  and  the  IBM  360  com- 
puter were  extensively  employed.   Application  of  these  power- 
ful computational  resources  to  the  available  data  would 
hopefully  statistically  reveal  what  parameters  are  most 
significant  in  marine-fog  prediction. 

B.  DATA  DESCRIPTION 

Upper-air  soundings  and  surface  3-hourly  airways  reports 
for  the  five-year  period  of  1  July  1968  to  30  June  1973  were 
used  for  three  U.  S.  west  coast  stations:   San  Diego,  Cali- 
fornia; Oakland,  California;  and  Quillayute,  Washington. 
This  initial  data  base  included  almost  11,000  soundings  (two 
daily)  and  over  43,700  surface  reports  (eight  daily).   Sea- 
surface  temperature  information  for  the  period  of  interest 
was  extracted  from  the  Fishing  Information  Bulletin  Supple- 
ments published  by  the  U.  S.  Department  of  Commerce.   The 
data  provided  by  the  Naval  Weather  Service  Detachment, 


23 


Asheville,  North  Carolina  were  taken  from  the  National  Cli- 
matic Center's  historical  data  files,  known  as  Tape  Data 
Family-14  and  -56  (TDF-14  and  TDF-56),  Airways  Surface  Obser- 
vations, and  Radiosonde  Data. 

C.    DEFINITION  OF  FOG  AND  FOG  DAY 

For  purposes  of  this  study  "fog"  (as  in  marine  advection 
fog)  is  defined  as  a  visible  aggregate  of  minute  particles 
of  water  (droplets)  based  at  the  earth's  surface,  which  re- 
duces horizontal  visibility  to  less  than  seven  mi les  and  is 
an  obstruction  to  vision  at  the  reporting  station,   A  "fog 
day"  as  used  in  this  study  was  defined  as  the  period  from 
1601  PST  on  one  calendar  day  through  1600  PST  on  the  follow- 
ing calendar  day,  so  that  a  day  would  usually  include  both 
the  formation  and  dissipation  time  of  fog  (after  Leipper, 
1948)  . 
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IV.    EXPERIMENTAL  PROCEDURES 

A.    DATA  PREPARATION  AND  ORGANIZATION 

In  any  study  involving  large  amounts  of  information  per- 
haps the  most  difficult  and  crucial  first  step  is  the  proper 
preparation  and  organization  of  the  data.   Initially  the 
data  tapes  had  to  be  screened  for  missing  reports  and  modi- 
fied for  the  use  of  FORTRAN  programming  on  the  NPS  IBM  360 
computer.   All  surface  and  upper-air  reports  were  identified 
by  year,  month,  day  and  hour  numeric  time  groups.   The  upper- 
air  sounding  time  group  identifiers  were  converted  from 
Greenwich  Mean  Time  (GMT)  to  Pacific  Standard  Time  (PST)  in 
order  to  conform  with  the  airways  surface  observation  time 
groups.   The  sea-surface  temperature  (SST)  data  were  extrac- 
ted twice  a  month  (1st  to  the  15th  and  16th  to  the  30th/31st) 
from  the  SST  charts  in  the  Fishing  Information  Bulletin  Sup- 
plements (U.  S.  Department  of  Commerce,  1968-1973),  and  were 
added  to  the  surface  report  data.   The  temperature  informa- 
tion on  the  surface  reports  was  then  converted  from  the 
Fahrenheit  scale  to  the  Celsius  scale,  thereby  utilizing  the 
metric  system  for  all  data  fields. 

Each  3-hourly  surface  observation  contains  35  data  fields 
(not  including  the  tape  and  station  number  and  the  date) , 
which  describe  the  following  meteorological  parameters: 
ceiling,  visibility,  wind  direction  and  speed,  dry-bulb 
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temperature,  wet-bulb  temperature,  dew-point  temperature, 
relative  humidity,  sea-level  pressure,  station  pressure,  sky 
condition,  cloud  information  (amount,  type,  and  height  of  as 
many  as  four  cloud  layers  plus  the  total  amount  and  total 
opaque  amount,  as  well  as  summation  amounts  at  the  second 
and  third  layers),  and  other  atmospheric  phenomena,  includ- 
ing wind  storms,  liquid  and  frozen  precipitation,  and  ob- 
structions to  vision.   This  last  data  field  is  especially 
important  since  the  definition  of  a  fog  occurrence  is  based 
upon  the  presence  of  fog  at  the  station  as  an  obstruction 
to  vision.   Furthermore,  the  code  lists  three  types  of  fog: 
"fog",  "ice  fog",  and  "ground  fog".   Although  marine  advec- 
tion  fog  is  not  explicitly  indicated,  the  great  majority  of 
fog  occurrences  were  described  as  "fog"  and  are  probably 
marine  fogs  considering  the  three  stations  selected  for  this 
study . 

The  upper-air  soundings  were  reported  twice  daily  (0400 
and  1600  PST)  and  contained  the  following  data  for  both  man- 
datory and  significant  levels:  pressure  in  millibars  (mb), 
height  in  meters  (m) ,  temperature  in  degrees  Celsius,  rela- 
tive humidity,  and  wind  direction  and  speed.   Since  fog  is 
a  surface  phenomenon,  the  entire  sounding  was  not  needed. 
As  reported  by  McClure  (1974),  the  important  inversion  in- 
formation relating  to  the  occurrence  of  fog  existed  at  850  mb 
and  below.  Consequently,  each  upper-air  report  was  arbi- 
trarily truncated  at  700  mb . 

The  surface  observations  and  the  upper-air  reports  were 
integrated  on  one  tape  by  numeric  time  groups  and  organized 
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by  location.   The  five  years  of  data  were  subdivided  into  two 
seasons,  designated  "dry"  and  "wet".   This  classification  was 
prompted  by  the  fact  that  west  coast  weather,  especially 
California,  has  two  dominant  seasons  characterized  by  a 
generally  dry  summer  period  and  a  wet  winter  period.   The 
so-called  "dry"  season  consisted  of  the  six  months  from  April 
through  September  while  the  remaining  six  months,  October 
through  March,  comprised  the  "wet"  season.   Initial  attempts 
to  study  the  data  by  months  instead  of  by  seasons  was  unsuc- 
cessful because  the  monthly  sample  size  of  fog  cases  was 
often  too  small.    (This  was  especially  true  for  the  summer 
months  at  Oakland,  California.) 

Nineteen  meteorological  parameters,  derived  from  both 
the  surface  and  the  upper-air  data,  were  selected  for  analy- 
sis.  Since  upper-air  observations  were  made  only  twice  a  day, 
the  nineteen  parameters  were  calculated  for  the  times  of  the 
upper-air  soundings,  namely  0400  PST  and  1600  PST.   Neverthe- 
less, each  3-hourly  surface  report  was  checked  for  the  occur- 
rence of  fog  and  such  information  was  retained  in  the 
parameter  "fog  strength",  FS .    (A  discussion  of  the  parameters 
is  found  in  section  B.) 

As  the  parameters  were  calculated,  the  reports  were  fur- 
ther classified  into  four  categories,  based  on  the  occurrence 
or  non-occurrence  of  fog  over  the  previous  three  days,  and 
whether  fog  was  reported  within  the  current  24-hour  period. 
If  fog  (no  fog)  occurred  between  1601  PST  on  one  calendar 
day  and  1600  PST  on  the  following  calendar  day,  that  period 
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was  classified  as  a  "fog"  day  ("no-fog"  day).   If  fog  was 
reported  on  any  surface  observation  in  the  72  hours  (3  days) 
prior  to  the  current  fog/no-fog  day,  that  period  was  termed 
"fog  history".   If,  on  the  other  hand,  no  fog  was  reported 
during  the  last  72  hours,  a  "no-fog  history"  period  had  oc- 
curred.  Using  these  definitions,  the  four  categories  shown 
in  Table  I  were  derived.   Generally  speaking,  categories  1 
and  2  depict  the  possible  onset  of  fog  since  they  have  no- 
fog  histories,  while  categories  3  and  4  depict  the  possible 
di  sper s  ion  of  fog  since  they  have  fog  histories. 

The  surface  observations  and  upper-air  data  were  now  in 
a  form  which  could  be  easily  analyzed.   Application  of  the 
discriminant  analysis  technique  required  the  data  to  appear 
in  two  or  more  groups.   The  categories  with  common  histories 
and  opposite  current  fog  designations  were  concatenated  and 
fed  into  BMD07M  in  order  to  develop  discriminant  functions. 
For  example,  categories  1  and  2  (3  and  4)  were  used  together 
where  group  1(3)  consists  of  "foggers"  and  group  2(4)  con- 
sists of  "no-foggers"  in  the  current  24-hour  period. 

1 ,   Data  Limitations  and  Special  Considerations 

In  any  set  of  data  there  are  always  certain  problems 
or  limitations  which  must  be  reconciled  or  at  least  recognized 
The  most  obvious  restriction  in  this  study  was  that  the  data 
covered  only  a  five-year  period  (July  1958  to  June  1973)  and 
is  given  for  just  three  west  coast  stations.  These  stations 
are  separated  by  hundreds  of  miles  and  extend  over  15  degrees 
of  latitude.   Also,  recall  that  the  surface  observations  were 
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made  at  three-hour  intervals  (0100,  0400,  0700,  1000,  1300, 
1600,  1900  and  2200  PST) ,  while  the  upper-air  soundings  were 
launched  only  twice  a  day  (0400  and  1600  PST).   It  is  impor- 
tant to  realize  that  although  fog  may  have  occurred  during 
a  24-hour  period,  say  at  1000  PST,  a  full  list  of  the  para- 
meters used  for  fog  discrimination  was  available  only  at  the 
sounding  times,  0400  and  1600  PST.   Considerable  amounts  of 
fog  were  reported  at  times  intermediate  to  the  upper-air 
data.   This  implies  that  significant  changes  in  the  existing 
sounding,  especially  in  the  marine  layer,  may  have  occurred 
in  the  interim,  resulting  in  the  formation  of  fog.   These 
subtle  changes  would  not  be  reflected  necessarily  in  either 
the  upper-air  or  surface  observations  for  that  fog  day. 

Several  other  limitations  of  the  data  set  should  be 
mentioned.   First  of  all,  if  any  information  field  in  either 
the  surface  observations  or  the  radiosonde  soundings  was 
missing  in  a  reporting  period  (0400  or  1600  PST),  that  re- 
port had  to  be  discarded.   Unfortunately,  missing  data  or  a 
missing  report  caused  a  gap  in  the  72  hours  of  "fog  history" 
and  resulted  in  a  loss  of  three  days  of  information.   Never- 
theless, only  about  five  percent  of  all  the  reports  were 
rejected  as  a  result  of  missing  data.   Another  restriction 
which  the  data  set  imposed  upon  this  research  was  that  the 
broad  synoptic  picture  was  not  considered.   That  is,  sounding 
stations  were  not  used  to  derive  advective  quantities  or 
define  air-mass  boundaries,  nor  were  synoptic  charts 
employed.   Furthermore,  the  sea-surface  temperature  (SST) 
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data  taken  from  the  Fishing  Information  Bulletin  Supplements 
may  not  always  be  the  most  representative  SST  nearest  the 
surface  station.   Although  the  SST  changes  rather  slowly,  a 
bi-weekly  average  SST  can  be  inaccurate,  especially  in  the 
summer  months  when  solar  radiation  can  result  in  significant 
temperature  changes.   The  Oakland  station  SST  was  a  special 
problem  because  the  station  is  located  on  the  east  side  of 
San  Francisco  Bay,  several  miles  from  the  ocean  regime.   The 
sea-surface  temperatures  nearest  the  bay  outlet  were  used, 
since  the  marine  fog  which  forms  offshore  and  affects  Oak- 
land generally  enters  through  this  opening. 

One  final  special  consideration  was  recognized  in 
San  Diego  where  the  surface  observation  station  (Lindberg 
Field)  and  the  upper-air  station  (Montgomery  Field)  are  not 
colocated.   The  radiosonde  launch  site  is  about  5.4  miles 
inland  from  Lindberg  Field  and  is  elevated  124  meters  above 
sea  level  (or  about  115  meters  above  the  surface  station). 
Under  these  circumstances  a  low-level  inversion  within  the 
124-meter  layer  could  exist  and  yet  not  be  shown  on  the 
sounding.   McClure  felt  that  in  fog  situations  where  inver- 
sions did  not  exist,  but  seemingly  should  have,  there  may 
have  been  an  inversion  near  the  sea  surface.   Consequently, 
in  this  study  a  special  procedure  was  implemented  to  "bogus 
in"  a  surface  inversion  under  the  proper  circumstances.   If 
the  surface  temperature  at  Montgomery  Field  (upper-air 
station)  was  equal  to  or  greater  than  the  surface  tempera- 
ture at  Lindberg  Field  (surface  observation  station)  an  iso- 
thermal condition  or  low-level  inversion  was  said  to  have 
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existed.   If  the  sounding  had  no  inversion  shown,  the  thick- 
ness of  the  low-level  inversion  was  taken  to  be  the  differ- 
ence in  elevation  between  Montgomery  and  Lindberg  Fields 
(124m-9m=115m) .   It  is  recognized  that  this  is  an  artificial- 
ity and  that  other  natural  influences  may  account  for  these 
differences,  such  as  the  nature  of  the  underlying  surface 
and  the  difference  between  the  marine  air  and  urban  environ- 
ments . 

B.    DISCUSSION  OF  THE  FOG/NO-FOG  DISCRIMINATION  PARAMETERS 

The  selection  of  the  meteorological  parameters  to  be 
used  as  fog  discrimination  variables  was  facilitated  by  pre- 
vious researchers  in  this  subject.   The  three  non-diurnal 
indices  developed  by  Leipper  (1948)  as  fog  forecasting  aids 
in  San  Diego  were  adopted  with  only  minor  modifications  for 
use  in  this  studyo   McClure's  thesis  (1974)  provided  informa- 
tion about  three  important  inversion  parameters:  height  of 
the  inversion  base  (also  one  of  Leipper's  indices),  thickness 
of  the  inversion  layer,  and  the  temperature  gradient  within 
the  inversion  layer. 

Using  the  works  of  Leipper  and  McClure  as  guidance,  nine- 
teen initially  selected  parameters  were  chosen  or  computed 
from  the  surface  observations  and  upper-air  data.   These 
parameters  are  defined  as  follows  (all  temperatures  in  degrees 
Celsius,  heights  in  meters,  wind  speeds  in  knots  and  pressures 
i  n  mi llibar s )  : 

1.   TB :   the  temperature  at  the  base  of  the  inversion  if 
an  inversion  exists,  otherwise  zero. 
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TT: 


SI: 


HB: 


the  temperature  at  the  top  of  the  inversion  if 
an  inversion  exists,  otherwise  zero. 
the  strength  of  the  inversion,  measured  as  the 
difference  between  the  temperature  at  the  top 
of  the  inversion  and  the  temperature  at  the 
base  of  the  inversion  (TT-TB=SI),  if  an  inver- 
sion exists,  otherwise  zero, 
the  height  of  the  base  of  the  inversion  if  an 


inversion  exists,  otherwise  zero. 

5.  HT:   the  height  of  the  top  of  the  inversion  if  an 

inversion  exists,  otherwise  zero. 

6.  THK:  the  thickness  of  the  inversion,  measured  as  the 

difference  between  the  height  of  the  top  of  the 
inversion  and  the  inversion  base  (HT-HB=THK), 
if  an  inversion  exists,  otherwise  zero. 

7.  DIR:  the  direction  from  which  the  wind  is  blowing 

at  the  950-mb  level. 

8.  SPD:  the  wind  speed  at  the  950-mb  level. 

9.  TI :   the  temperature  index,  measured  as  the  differ- 

ence between  the  temperature  at  the  top  of  the 
inversion  (or  the  surface  dry-bulb  temperature 
if  there  is  no  inversion)  and  the  sea-surface 
temperature  (TT(or  TDB) -SST=TI ) „ 
10.   WDB :  the  surface  wet-bulb  depression,  measured  as 

the  difference  between  the  surface  dry-bulb  and 
wet-bulb  temperatures  (TDB-TWB=WBD) . 
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11.  DPD:  the  dew-point  depression,  measured  as  the  dif- 

ference between  the  surface  dry-bulb  and  dew- 
point  temperatures  (TDB-TDP=DPD) . 

12.  WET:  the  surface  moisture  index,  measured  as  the 

difference  between  the  dew-point  and  the  sea- 
surface  temperatures  ( TDP-SST  =  WET)  . 

13.  TDB :  the  surface  dry-bulb  temperature. 

14.  TWB :  the  surface  wet-bulb  temperature. 

15.  TDP :  the  surface  dew-point  temperature. 

16.  SST:  the  sea-surface  temperature. 

17.  RH :   the  surface  relative  humidity. 

18.  SLP:  the  sea-level  pressure. 

19.  FS :   the  fog  strength;  an  integer  number  from  0  to 

4  indicating  the  number  of  three-hourly  surface 
observations  during  a  twelve-hour  period  (from 
0400-1600  PST  or  1600-0400  PST)  reporting  fog. 
The  notation  used  to  express  the  parameters  and  their 
relative  day  and  time  within  the  three-day  history  period  is 
quite  simple.   The  pre- subs cript  denotes  the  day  during  this 
period,  while  the  post-subscript  denotes  the  hour  of  observa- 
tion (0400  or  1600  PST).   Examples  of  this  notation  are  as 

f ol lows : 
H  T 

1   16  =  the  height  of  the  top  of  the  inversion,  taken 
from  the  first  day  of  history  and  the  1600  PST 
report.    (This  would  be  the  most  recent  observa- 
tion in  the  history  period.) 
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SI 
3   04  =  the  strength  of  the  inversion,  taken  from  the 

third  day  of  history  and  the  0400  PST  report. 

(This  would  be  the  oldest  observation  in  the 

history  period.) 

In  addition  to  the  nineteen  basic  parameters  described 

above,  12-,  24-,  and  48-hour  time  di  f  f er ences  in  each  of 

these  variables  were  calculated  using  the  three  days  of 

"history"  associated  with  a  fog  occurrence.   Because  of  the 

BMD07M  program  restriction  of  80  variables  there  could  be 

computed  only  76  (of  a  possible  95)  12-hour  differences,  76 

24-hour  differences  and  38  48-hour  differences  from  the 

total  list  of  114  (i.e.,  6x19)  parameters  for  each  day. 

These  differences  can  be  shown  schematically  as  follows: 

12-hour  differences 

-1600  (19  variables) 
=  0400 

=16  0  0  " . 

=0400  " 

=1600  " 

-0400  " 


3  days  of 
history  as- 
sociated 
with  a  fog 
occur renc  e 


5x19=95  possible  differences 

24-hour  differences 
rl600   (19  variables) 
0  4  0  0         " 


2:  1=1600 
0400 

4:  ^1600 
0400 


4x19=76  possible  differences 
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48-hour  differences 
pieOO   (19  variables) 
0400 


1600 

0400 
1600 
0400 


^2'^^6 


2x19=38  possible  differences 
Examples  of  the  notation  used  to  represent  these  differences 
are  : 

A  TB,^  ^^  =  the  first  (1600-0400)  12-hour  difference 

1    1 D- 04 

in  the  temperature  at  the  base  of  the  in- 
vers  ion . 

the  second  (1600-1600)  24-hour  difference 
in  the  temperature  at  the  base  of  the  in- 
version. 
A  TB    =  the  (0400-0400)  48-hour  difference  in  the 
temperature  at  the  base  of  the  inversion. 
The  time  differences  in  the  nineteen  basic  parameters 
are  considered  to  be  important  since  they  allow  for  the  in- 
fluence of  persistence,  advection  and  other  physical  changes 
in  the  variables.   Thus,  the  analyst  can,  at  least,  theoreti- 
cally, project  "trends"  (time  differences)  of  the  parameters 
ahead  in  time  for  the  purpose  of  developing  a  forecast. 

C.    DISCRIMINANT  ANALYSIS 

The  use  of  linear  discriminant  analysis  in  this  study  has 
been  of  fundamental  importance  in  the  successful  development 


Taken  from  Chapter  23  (written  by  VJilliam  R.  Klecka)  of 
Statistical  Package  for  the  Social  Sciences,  Norman  H.  Nie, 
et  al,  McGraw-Hill,  1970. 
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of  classification  equations  for  marine  fog  prediction.   The 
objective  of  the  analysis  is  to  statistically  distinguish 
between  two  or  more  categories  of  cases.   In  this  research 
these  categories  (see  Table  I)  were  determined  by  the  occur- 
rence or  non-occurrence  of  fog  during  three  successive  days, 
and  whether  fog  was  reported  during  the  current  24-hour 
period.   Thus,  the  groups  "fog"  and  "no-fog"  were  established 
as  previously  described. 

In  order  to  distinguish  between  the  fog  and  no-fog 
groups,  the  nineteen  initially  selected  parameters  and  their 
time  differences  were  used  as  "discriminating  variables"  to 
measure  characteristics  on  which  these  groups  are  expected 
to  differ.   The  mathematical  objective  of  discriminant  analy- 
sis is  to  weight  and  linearly  combine  the  discriminating 
variables  in  some  manner  so  that  the  groups  are  forced  to 
be  as  statistically  distinct  as  possible.   In  other  words, 
one  hopes  to  be  able  to  "discriminate"  between  the  groups 
in  the  sense  of  being  able  to  tell  them  apart. 

Discriminant  analysis  attempts  to  accomplish  this  by 
forming  one  or  more  linear  combinations  of  the  discriminat- 
ing variables.   These  "discriminant  functions"  are  of  the 
form : 


D   =  c  ^  +  c  ,V^  +  c  ^V^  +  ...  +  c   V 

m     mO     ml  1     m2  2  mn  n 


(1) 


where   D    is  the  discriminant  score  for  group  m,  the 


m 


mj 


are  the  coefficients  with   c  ^   being  the  constant,  and  the 

mO 

V's   are  the  raw  values  of  the   n   discriminating  variables 
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used  in  the  analysis.  There  is  always  a  separate  equation 
for  each  group.  Hopefully,  the  discriminant  scores  (D's) 
for  the  cases  within  a  particular  group  will  be  fairly  sim- 
ilar. Nevertheless,  the  functions  are  formed  in  such  a  way 
as  to  maximize  the  separation  of  the  groups.  Once  the  dis- 
criminant functions  have  been  derived,  the  primary  research 
objective  of  classification  can  be  pursued. 

Classification  is  the  process  of  identifying  the  group 
to  which  a  fog/no-fog  case  is  most  likely  to  belong  when  the 
only  information  known  is  the  case's  values  on  the  discrim- 
inating variables.   The  classification  procedure  involves 
the  use  of  a  separate  linear  combination  of  the  discriminat- 
ing variables  for  each  group.   These  produce  a  probability 
of  membership  in  the  respective  group,  and  the  case  is  as- 
signed to  the  group  with  the  highest  probability  of  occur- 
rence.  One  problem  with  this  type  of  classification  proce- 
dure is  that  the  rule  of  highest  probability  defines  a  very 
strict  dividing  line.   A  .51  probability  of  a  fog  occurrence 
versus  a  .49  probability  of  no  fog  would  necessarily  lead 
to  a  fog  classification,  yet  the  situation  is  not  really  all 
that  clear.   A  Bayesian  adjustment  of  the  posterior  proba- 
bilities is  often  desirable  when  the  costs  of  misclassifica- 
tion  into  certain  groups  are  high,  or  when  the  groups  are  of 
grossly  different  sizes.   In  this  study  the  posterior  proba- 
bilities were  not  adjusted. 

During  the  analysis  phase  of  this  study  the  primary  use 
of  classification  was  to  test  the  adequacy  of  the  derived 
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discriminant  functions.   By  classifying  the  fog  and  no-fog 
cases  used  to  derive  the  functions  in  the  first  place  and 
comparing  predicted  group  membership  with  actual  group  mem- 
bership, one  can  empirically  measure  the  success  in 
discrimination  by  observing  the  proportion  of  correct 
classifications . 

The  BMD07M  program  used  in  this  research  is  a  s  tepwi  se 
linear  discriminant  analysis  routine  which  can  accept  a 
maximum  of  80  discriminating  variables.   Since  there  are 
generally  more  variables  available  than  are  necessary,  the 
stepwise  analysis  procedure  is  very  helpful.   It  begins  by 
choosing  the  single  best-discriminating  variable  as  deter- 
mined by  the  selection  criterion.   A  second  discriminating 
variable  is  selected  as  the  variable  best  able  to  improve 
the  value  of  the  discrimination  in  combination  with  the 
first.   The  third  and  subsequent  variables  are  chosen  in  a 
similar  manner  according  to  their  ability  to  contribute  to 
further  discrimination.   At  each  step,  variables  already 
selected  may  be  removed  if  they  are  found  to  reduce  discri- 
mination when  combined  with  the  more  recently  selected 
variables.   This  process  continues  until  the  remaining  vari- 
ables are  no  longer  able  to  contribute  significantly  to 
further  discrimination. 

The  stepwise  analysis  procedure  just  described  has  sev- 
eral test  criteria  which  control  the  entry  or  removal  of 
variables  from  the  discrimination  equationo   Ultimately, 
the  selection  criterion  is  determined  by  the  F  ratio  or 
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F  value .   The  F  value  is  the  likelihood  ratio  of  equality  on 
the  test  variable  over  all  the  groups  (two  groups  in  this 
study),  given  the  distribution  produced  by  the  variables  al- 
ready entered.   In  other  words,  the  "F-to-enter"  value  is  a 
test  for  the  statistical  significance  of  the  amount  of  cen- 
troid  separation  added  by  this  variable  above  and  beyond  the 
separation  produced  by  the  previously  entered  variables.   At 
each  step  in  the  BMD07M  program  there  is  a  similar  test  made 
of  all  variables  already  selected  which  have  an  inclusion 
level  of  one  (the  inclusion-level  default  option  of  one  was 
used  throughout  this  study) .   Here,  the  test  is  whether  the 
particular  variable  still  adds  a  significant  amount  to  the 
group  separation,  given  the  other  variables  now  in  the  equa- 
tion.  As  more  variables  are  chosen,  it  is  possible  that 
some  of  those  entered  earlier  will  no  longer  be  contributing. 
This  occurs  because  the  information  that  they  contain  about 
group  differences  is  now  available  in  some  combination  of 
the  other  included  variables.   The  contribution  is  measured 
by  the  "F- to-r emove " ,  which  is  an  F  test  of  the  discrimina- 
tion power  currently  added  by  the  variable  in  question.   If 
a  variable  already  entered  has  an  F-to-remove  value  less  than 
one  (again,  the  default  option  was  used),  it  will  be  removed 
from  the  equation. 

The  BMD07M  program  computes  unstandardized  coefficients 
in  the  classification  functions,  thus  making  analysis  and 
interpretation  of  such  coefficients  difficult.   The  unstand- 
ardized coefficients  do  not  report  the  relative  importance 
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of  the  variables  since  they  have  not  been  adjusted  for  the 
measurement  scales  and  variability  in  the  original  variables. 
However,  the  uns tandardized  discriminant  function  coeffi- 
cients can  be  easily  utilized  to  find  the  fog/no-fog  proba- 
bility.  In  order  to  do  so  one  first  computes  the 
discriminant  score  by  multiplying  the  coefficients  by  the 
raw  values  of  the  associated  variables,  summing  them  together, 
and  adding  a  constant  to  adjust  for  the  grand  means.   This 
process  is  represented  by  Eq.   (1) .   Then  the  posterior  proba- 
bility of  case  k  having  come  from  group  m  is: 


p   exp(D    ) 
m m  ,  K 


m,k    p^  ^^P(Df^,^)  +  P„f  exp(D^^^j^) 


(2) 


where    m  =  a  fog  or  no-fog  day  (f  or  nf) 

p  =  the  prior  probability  of  group  m 
m 

D   , =  the  discriminant  score  for  case  k  of  group  m. 
m ,  k 

In  this  study  the  prior  probabilities  of  belonging  to  either 
the  fog  or  no-fog  group  were  taken  to  be  equal,  so  p   was 
0.5.   Some  discriminant  analysis  programs  allow  for  automa- 
tic adjustments  to  the  prior  probabilities  on  the  basis  that 
such  probabilities  are  proportional  to  the  number  of  cases 
in  each  group,  i.e.,  cases  will  be  more  likely  assigned  to 
a  larger  group.   The  BMD07M  routine  does  not  have  an  auto- 
matic adjustment  feature,  but  manual  changes  in  the  prior 
probabilities  can  be  accomplished. 

A  wealth  of  statistical  information  can  be  derived  from 
the  BMD07M  computer  program--much  more  than  this  researcher 
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required.   The  information  generally  used  in  this  study  con- 
sisted of  : 

1.  Group  means  and  standard  deviations 

2 .  At  each  step : 

a.  Variables  included  and  F  value  to  remove 

b.  Variables  not  included  and  F  value  to  enter 

3.  After  each  of  the  first  seven  steps  and  after  the 
last  step  the  classification  functions  were  printed 

4.  The  posterior  probability  of  coming  from  each  group 
was  computed  for  each  case 

5.  Summary  table.   For  each  step  of  the  procedure  the 
following  were  utilized: 

a.  Variable  entered  or  removed 

b.  F  value  to  enter  or  remove 

c.  Number  of  variables  included 

D.   ANALYSIS  ROUTINE 

After  the  data  had  been  organized  by  location,  season 
and  fog  category,  the  job  of  analysis  was  finally  initiated. 
Since  considerable  fog  research  has  been  done  in  the  San 
Diego  area  (see  Section  II),  this  study  began  here. 

For  each  location  and  for  both  the  "dry"  and  "wet"  sea- 
sons, the  fog  and  no-fog  cases  of  categories  1  and  2  (and 
later  3  and  4)  were  combined  and  read  into  the  BMD07M  rou- 
tine.  It  is  important  to  remember  the  program  restriction 
of  a  maximum  of  80  discriminating  variables,  since  this  lim- 
itation influenced  the  analysis  methodology.   Associated  with 
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each  fog  or  no-fog  occurrence  were  3  days  (six  reporting 
periods)  of  history,  listing  a  total  of  114  variables  (6 
reports  x  19  parameters  =  114  variables).   Also  recall  that 
there  were  the  12-,  24-,  and  48-hour  time  differences  of 
these  variables  which  were  available,  thus  adding  another 
209  possible  fog  discriminators.   However,  for  any  one  com- 
puter program  no  more  than  80  of  these  variables  could  be 
utili  zed . 

In  order  to  test  the  discriminating  potential  of  all  the 
possible  variables,  seven  separate  analysis  programs  were 
completed  for  each  location,  season  and  category  pair  (1  and 
2;  3  and  4) .   These  programs  differed  from  one  another  by 
virtue  of  the  variables  being  read  into  the  BMD07M  routine 
and  are  described  as  follows: 

1.  The  19  initially  selected  parameters  for  only  the 
first  two  days  of  history,  but  both  1600  and  0400 
PST  reports  were  used. 

2.  The  19  basic  parameters  for  all  3  days  of  history, 
but  only  the  1600  PST  reports  were  used. 

3.  The  19  basic  parameters  for  all  3  days  of  history, 
but  only  the  0400  PST  reports  were  used. 

4.  Only  the  12-hour  time  differences  of  the  19  basic 
parameters  were  used  (four  12-hour  differences). 

5.  The  24-hour  time  differences  for  only  the  1600  PST 
reports  were  used. 

6.  The  24-hour  time  differences  for  only  the  0400  PST 
reports  were  used. 
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7,   The  48-hour  time  differences  for  both  reporting 
periods  were  used. 

Each  of  these  computer  programs  executed  the  discriminant 
analysis  routine  with  the  available  variables  and  printed 
a  summary  table  listing  the  discriminating  variables  in 
order  of  entrance  and  the  classification  functions.   After 
completing  these  jobs  for  both  seasons  (and  for  a  given  fog 
category  pair),  the  best  seven  or  eight  discriminating 
variables  were  selected  from  each  program  and  put  in  a  list. 
Any  of  these  variables  which  may  have  been  duplicated  be- 
tween, say,  programs  #1  and  #3  (see  list  above) ,  or  between 
wet  and  dry  seasons  were  discarded.   Finally,  the  "best"  80 
discriminating  variables  were  computed  in  a  "composite" 
program  for  final  analysis.   This  composite  program  was 
executed  for  both  seasons  and  resulted  in  the  final  classi- 
fication functions  for  a  given  location  and  fog  category 
pair.   Then  the  entire  process  was  repeated  for  the  other 
fog  category  pair  (3  and  4) .   Thus,  for  each  location  a 
minimum  of  32  discriminant  analysis  programs  v;ere  executed 
in  order  to  build  the  final  classification  equations.   When 
using  these  equations  one  can  calculate  the  percentage  pro- 
bability of  fog/no  fog  over  a  24-hour  period  beginning  at 
1601  PST  on  the  current  calendar  day. 
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V.   ANALYSIS  RESULTS 

A.    INTRODUCTORY  REMARKS 

The  analysis  results  of  this  research  are  presented 
largely  in  tabular  fornio   For  each  west  coast  station  studied 
there  are  three  different  types  of  tableSo   The  first  type 
shows  the  number  of  usable  fog  and  no-fog  "cases"  (sometimes 
called  "foggers"  or  "no-f oggers " )  for  each  analysis  category, 
A  usable  fog  (no-fog)  "case"  is  simply  a  fog  (no-fog)  day 
for  which  there  are  a  complete  list  of  parameterso   The  second 
type  of  table  presents  the  classification  functions,  the 
variables  utilized,  and  their  associated  F  values.   Recall 
that  the  F  value  (F  ratio)  shows  the  relative  importance  of 
the  variable  in  question  to  the  existing  equationo   The  F 
values  shown  in  the  tables  are  those  for  the  final  classifi- 
cation equation„   Thus,  in  general,  the  larger  the  F  value 
(relative  to  unity)  for  a  given  variable,  the  greater  that 
variable  contributes  to  the  discrimination  (centroid  separa- 
tion) between  the  fog  and  no-fog  groups.   The  third  type  of 
table  shows  the  season  and  category  for  a  given  location, 
and  the  percentage  of  fog  and  no-fog  cases  correctly  clas- 
sified by  the  discriminant  functions. 

It  is  important  to  remember  that  the  stepwise  dis- 
criminant analysis  procedure  results  in  an  optimal  set  of 
variables  being  selected  and  that  the  initially-best  dis- 
criminating variables  are  chosen  first.   Therefore,  in  the 
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tables  presenting  the  classification  functions  the  variables 
are  listed  in  order  of  decreasing  importance  according  to 
initial  entry.   Where  a  variable  had  been  entered  and  was 
subsequently  removed  there  appears  a  remark  indicating  the 
step  in  the  program  that  the  variable  was  deleted. 

B.    SAN  DIEGO  RESULTS 

Tables  II  and  III  show  the  number  of  fog  and  no-fog 
cases  for  the  San  Diego  dry  and  wet  seasons,  respectively. 
During  the  dry  season  the  fog  cases  comprise  only  15.5%  of 
the  total  number  of  cases,  while  during  the  wet  season  over 
28%  were  fog  days. 

Table  IV  presents  the  fog  and  no-fog  classification  func- 
tions for  the  dry  season,  categories  1  and  2.   The  F  values 
indicate  that  the  two  most  important  discriminating  variables 
in  the  final  equations  are  the  temperature  index  (hTI^^.)  and 

the  inversion  thickness  (,THK^^)  taken  from  the  most  recent 
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0400  PST  observations.   In  this  case  the  mean  . TI    for  the 
foggers  was  approximately  4.0  C,  while  the  no-foggers  had  a 
mean  of  about  1.6  C.   This  implies  that  a  relatively  large 
temperature  difference  between  the  air  temperature  (at  the 
surface  or  inversion  top)  and  the  SST  prevailed  just  prior  to 
most  fog  occurrences.   This  result  generally  supports  Leipper's 
(1948)  temperature  index  range  values  found  in  San  Diego. 
The  mean  values  of  the  inversion  thickness  for  the  fog  and 
no-fog  cases  differed  only  slightly  (385.7  m  versus  362.3  m), 
but  the  thickness  was  greater  for  the  fog  reports. 
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The  classification  functions  for  the  wet  season,  catego- 
ries 1  and  2,    are  given  in  Table  Vo   Although  the  moisture 
index  (  WET   )  and  the  inversion  strength  (^SI   )  were  en- 
tered as  the  first  two  variables,  they  were  subsequently  re- 
moved.  When  they  were  later  re-entered  they  made  only  a 
small  contribution  to  the  overall  discriminating  power  of 
the  final  equations.   In  this  program  the  ultimate  best  three 

variables  are  the  wet-bulb  depression  (.WBD.^),  the  height  of 
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the  top  of  the  inversion  (  HT   ) ,  and  the  temperature  index 
{  TI   )  at  1600  PST.   The  mean  wet-bulb  depression  for  the 
no-foggers  was  almost  lo3  degrees  greater  than  that  of  the 
fog  reportSo   The  three-day-old  HT  value  was,  on  the  average, 
81  m  higher  for  the  no-fog  caseSo   The  TI  mean  values  dif- 
fered by  about  1.1  C  (3o3  versus  2.2),  and  again  this  value 
was  usually  larger  for  the  fog  reports. 

Table  VI  shows  the  classification  functions  for  San  Diego 
during  the  dry  season  for  categories  3  and  4.   Recall  that 
categories  3  and  4  differ  from  1  and  2  in  that  they  have  a 
history  of  fogo   Furthermore,  categories  3  and  4  should  give 
some  clue  as  to  the  parameters  of  interest  in  fog  dispersion. 
The  F  values  once  more  identify  the  temperature  index  ^■\'^^rsA) 
as  being  very  important  to  the  final  equations,  and  surpris- 
ingly, so  is  the  three-day-old  wind  direction  (  DIR   ) o   As 
was  seen  before,  the  mean  TI  was  significantly  higher  for 
the  foggers  than  the  no-foggerso   In  this  instance  the  dif- 
ference was  about  2.5  degrees  (i.e.,  5o5  C  for  fog  versus 
3.0  C  for  no-fog) o   In  San  Diego  the  temperature  index   ap- 
pears to  be  an  important  fog/no-fog  discrimination  parameter. 
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The  950-mb  wind  direction  mean  values  indicate  that 
three  days  before  a  fog  occurrence  the  wind  was  southeasterly 
(about  126°),  whereas  prior  to  no-fog  days  the  wind  was  gen- 
erally from  the  southwest  (207°).   It  should  be  noted,  how- 
ever, that  the  standard  deviation  for  these  wind  directions 
was  about  120  degrees.   Therefore,  solid  conclusions  as  to 
the  significance  of  this  result  cannot  be  made  without  further 
analysis  on  a  broader  data  base.   Other  important  discriminat- 
ing parameters  shown  in  Table  VI  are  a  48-hour  change  in  the 

1600  PST  wet-bulb  temperature  (ATWB,^),  a  48-hour  change  in 
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the  0400  PST  wet-bulb  depression  (AWBD   ) ,  and  a  12-hour 

change  in  the  temperature  index  (A^TI^  .  ^  _)  . 

1   0  4-15 

The  San  Diego  wet  season  classification  functions  for 

categories  3  and  4  are  presented  in  Table  VII.   The  most 

capable  discriminators  in  these  equations  are  the  wet-bulb 

depression  (  WBD   ) ,  the  inversion-base  height  (  HE   ) ,  and 

the  fog  strength  (  FS,^),  all  measured  at  1600  PST.   The 
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statistics  for  the  wet-bulb  depression  mean  values  are  nearly 
identical  with  those  mentioned  for  categories  1  and  2  during 
the  wet  season.   The  no-fog  reports  averaged  over  1.3  degrees 
higher  than  those  of  the  fog  cases  (4.3  C  versus  3.0  C  res- 
pectively).  This  result  appears  to  be  consistent  for  all 
categories  in  the  San  Diego  winter  season.   The  mean  values 
of  the  inversion-base  height  exhibit  dramatic  differences  be- 
tween the  fog  and  no-fog  cases.   In  fog  situations  the  base 
heights  were  about  168  m,  whereas  the  heights  for  no-fog 
situations  were  535  m.   Thus,  as  the  fog  regime  begins  to 
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dissipate,  the  base  of  the  inversion  rises.   This  is  in 
agreement  with  Leipper's  model  for  San  Diego  winter  fogs 
(Leipper ,  1948)  . 

The  fog  strength  (FS)  parameter  can  be  a  significant 
discriminator  only  in  categories  3  and  4  where  there  is  a 
history  of  fog.   In  this  case,  during  the  wet  season  the 
mean  FS  value  for  1600  PST  was  nearly  0.88  for  fog  cases 
and  0.26  for  no-fog  reports.   This  simply  suggests  that  fog 
persi  s  tence  is  likely  to  occur  in  winter.   That  is,  fog  is 
likely  to  occur  in  one  of  the  three-hourly  surface  observa- 
tions in  successive  twelve-hour  periods. 

Table  VIII  shows  the  San  Diego  fog  and  no-fog  discrim- 
ination percentages  for  the  seasons  and  category  pairs.   The 
discriminating  capability  of  the  category  3  and  4  equations 
is  better  than  those  for  1  and  2,  especially  for  fog  classi- 
fication.  Overall,  the  weighted  average  fog  discrimination 
percentage  was  about  82.2,  while  the  no-fog  percentage  was 
nearly  81.3. 

C.    OAKLAND  RESULTS 

The  occurrence  of  fog  in  Oakland's  summer  (dry)  season 
is  indeed  infrequent.   For  the  five-year  data  period  of  this 
study  there  were  only  a  total  of  32  fog  days  during  this 
season,  as  shown  in  Table  IX.   This  amounts  to  less  than 
four  percent  of  the  total  cases.   Fog  occurrences  in  the  wet 
season  (shown  in  Table  X)  were  more  common--about  24%  of  the 
cases  reported  fog. 
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Table  XI  presents  the  fog  and  no-fog  classification  func- 
tions for  the  Oakland  dry  season,  categories  1  and  2.   As  the 
F  values  reveal,  the  best-discriminating  variables  in  the 
final  equations  are   the  24-hour   change  in  the  1600  PST  in- 
version thickness  (A  THK^^),  a  similar  time  difference  in  the 
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relative  humidity  (A  RH   ),  and  the  second  (0400-1600)  24- 
hour  change  in  the  temperature  index  (A  TI   ) .   It  is  in- 
teresting to  note  that  the  most  important  variables  in  this 
equation  are  all  24-hour  time  differences .   With  the  exception 

of  one  wind  speed  variable  (  SPD   ) ,  the  nineteen  basic  para- 

Z.  16 

meters  themselves  contribute  very  little  to  the  fog/no-fog 
discrimination . 

The  24-hour  thickness  change  mean  values  showed  that 
there  usually  was  a  decrease  in  the  1600  PST  inversion  thick- 
ness just  prior  to  a  fog  occurrence.   The  A  THK,_  average 
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for  the  fog  reports  was  about  -132  m,  whereas  the  no-fog 
value  was  2.6  m  (a  negligible  increase).   Although  the  stan- 
dard deviations  for  the  two  groups  average  about  265  m,  the 
means  clearly  indicate  a  reduction  of  the  inversion  thickness 
which  was  an  unexpected  result.   Since  this  season  contained 
only  25  fog  reports,  analysis  of  a  larger  data  set  seems 
necessary  in  order  to  be  certain  of  this  thickness  relation- 
ship . 

Another  surprising  result  in  this  analysis  was  that  the 
24-hour  relative  humidity  change  (A  RH   )  for  the  fog  reports 
was  an  average  -0.18%.   Although  the  mean  change  is  very 
small,  one  would  normally  expect  a  moderate  increase  in  the 
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relative  humidity  prior  to  fog.   A  RH ^ ^  for  the  no-fog 
cases  was  almost  zero,  and  the  standard  deviation  for  both 
groups  exceeded  2.0, 

The  mean  values  of  the  24-hour  change  in  the  tempera- 
ture index  (A  TI   )  reflected  an  increase  (almost  1.5  C) 
for  the  fog  cases  and  a  very  slight  decrease  (-0.06  C)  for 
the  no-fog  cases.   This  result  appears  to  be  consistent  v/ith 
San  Diego  observations;  several  days  prior  to  fog  develop- 
ment the  inversion  intensifies  which  generally  leads  to  an 
increase  in  the  temperature  index. 

The  classification  functions  for  Oakland's  wet  season, 
categories  1  and  2,  are  given  in  Table  XII.   Here  the  first 
two  variables  entered  in  the  stepwise  discriminant  analysis 
remained  the  best  variables  in  the  final  equations.   Both 
of  these  parameters  v;ere  from  the  most  recent  1500  PST  ob- 
servation and  are  the  temperature  at  the  top  of  the  inver- 
sion (,TT,^),  and  the  relative  humidity  (.RH^^).   The  means 
116  116 

of  the  TT  values  were  8.7  C  for  foggers  and  5.5  C  for  the 
no-foggers.   From  this  result  one  might  infer  that  the  in- 
version is  either  thicker  or  more  intense,  or  possibly  both. 
The  RH  mean  value  for  the  fog  cases  was  almost  ten  percent 
higher  than  for  the  no-fog  cases  (73.5%  versus  63.5%). 

Of  tertiary  importance  as  a  discriminating  variable  is 
the  first  (0400-1600)  12-hour  change  in  the  temperature  of 

the  inversion  base  (A  TB..  ,,).   The  mean  values  for  the  fog 
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and  no-fog  cases  were  approximately  3.2  C  and  1.4  C,  respec- 
tively.  Although  the  inversion  base  showed  a  12-hour  warming 
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tendency  for  both  groups,  the  average  increase  was  more  than 
twice  as  great  for  the  foggers.   The  reason  for  this  warming 
tendency  at  the  base  of  the  inversion  is  not  fully  under- 
stood by  this  researcher,  especially  since  this  was  over  a 
12-hour  period  (0400-1600)  during  which  diurnal  cooling  us- 
ual ly  occurs . 

Table  XIII  shows  the  classification  functions  for  Oak- 
land during  the  dry  season  for  categories  3  and  4.  As  the 
F  values  indicate,  the  three  most  capable  discriminators  in 

these  equations  are  the  fog  strength  (^FS^,),  the  950-mb 
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wind  speed  (  SPD   ) ,  and  the  wet-bulb  temperature  (  TWB   ) . 
Curiously,  all  three  of  these  variables  were  taken  from  the 
1600  PST  reports  which  were  two  days  oldo   It  is  important 
to  note  that  there  were  only  seven  fog  reports  for  this  sea- 
son and  category,  compared  to  over  eleven  times  as  many  (78) 
no-fog  cases  . 

Interestingly,  the  fog  strength  mean  for  the  fog  cases 
(0.14)  was  less  than  half  of  the  mean  for  the  no-fog  cases 
(0.36)  .   This  is  just  the  opposi  te  that  one  might  expect  if 
fog  persistence  had  occurred.   Because  the  fog  and  no-fog 
cases  were  so  disproportionate,  perhaps  some  adjustment  in 
the  prior  probabilities  reflecting  this  imbalance  would 
moderate  or  change  this  result  and  remove  a  potential  sta- 
tistical bias  . 

The  wind  speed  parameter  showed  less  than  one  knot  dif- 
ference between  the  mean  values  of  each  group — 3.9  knots  for 
foggers  and  4.8  knots  for  no-foggers.   The  wet-bulb 
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temperature  averages  were  about  16.0  C  and  14.3  C  for  the  fog 
and  no-fog  cases,  respectively.   Once  again,  a  larger  data 
sample  of  reports  is  required  to  prove  whether  these  differ- 
ences are  really  significant. 

The  Oakland  wet  season  classification  functions  for 
categories  3  and  4  are  presented  in  Table  XIV.   The  best  dis- 
criminator in  these  equations  is  the  temperature  at  the  inver- 
sion top,  measured  at  1500  PST  (-TT,^)--the  same  as  for  the 
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category  1  and  2  analysis  during  this  season.   The  mean 
values  for  this  parameter   indicate  nearly  a  4.3  degree 
spread  between  the  groups — about  10.0  C  for  the  foggers  and 
5.7  C  for  the  no-f oggers .   The  identical  inferences  may  be 
made  since  the  tendency  was  equal  to  that  of  the  category  1 
and  2  result. 

Two  other  important  discriminators  in  this  program  are 
the  48-hour  change  in  the  1600  PST  950-mb  wind  speed  (ASPD   ), 
and  the  0400  PST  height  of  the  inversion  top  (  HT   ) .   The 
average  wind  speed  change  was  only  a  one-third  knot  decrease 
for  the  fog  cases  and  a  0.6  knot  increase  for  the  no-fog 
reports.   Nevertheless,  in  this  data  set  the  nearly  one  knot 
difference  was  found  to  be  statistically  significant  by  the 
BMD07M  program.   The  mean  values  for  the  height  of  the  in- 
version top  for  foggers  and  no-foggers  were  566  m  and  537  m, 
respectively.   Although  the  difference  between  the  group 
means  was  only  71  m,  the  tendency  for  a  lower  near-surface 
inversion  prior  to  fog  occurrences  is  evident. 
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Table  XV  shows  the  Oakland  fog  and  no-fog  discrimination 
percentages  for  the  seasons  and  category  pairs.   The  high 
discrimination  score  for  the  dry  season,  categories  3  and  4, 
was  largely  due  to  the  unusually  small  number  of  fog  and  no- 
fog  cases.   If  the  BMD07M  program  had  been  allowed  to  con- 
tinue until  the  F-to-enter  value  dropped  below  1.0  (as  was 
done  with  all  other  equations) ,  the  resultant  equation  would 
have  had  30  variables  (after  48  program  steps!)  and  could 
achieve  a  "perfect"  discrimination  score.   However,  thirty 
variables  was  considered  to  be  an  excessive  number  for  this 
small  data  set  and  the  program  was  limited  to  twenty-one 
steps  which  included  the  final  fifteen  variables  shown  in 
Table  XIII.   The  overall  weighted  fog  discrimination  percent- 
age was  above  77.4,  while  the  no-fog  percentage  was  76.8. 

D.    QUILLAYUTE  RESULTS 

Whereas  the  occurrence  of  fog  in  Oakland  was  a  rarity  in 
summer,  Quillayute,  Washington  has  fog  as  a  fairly  frequent 
visitor  throughout  the  entire  year.   In  fact,  for  the  period 
of  this  research  exactly  50%  of  the  dry-season  reports  and 
nearly  54%  of  the  wet-season  reports  were  fog  cases.   The 
number  of  fog  and  no-fog  cases  for  these  two  seasons  are 
shown  in  Tables  XVI  and  XVII.   It  is  likely  that  in  this  area 
a  significant  percentage  of  fog  occurrences  were  associated 
with  frontal  activity,  since  most  mid-latitude  storms  track 
across  the  Washington  coast. 
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Table  XVIII  presents  the  fog  and  no-fog  classification 
functions  for  the  Quillayute  dry  season,  categories  1  and  2, 
There  are  many  variables  with  relatively  large  F  values,  but 
the  three  best  discriminators  in  these  equations  are  the 
most  recent  24-hour  change  in  the  0400  PST  temperature  index 

(A  TI   ) ,  the  second  12-hour  change  (0400-1600  PST)  in  the 
temperature  index  (A  TI       ),  and  the  950-mb  wind  direction 
measured  at  0400  PST  (  DIR   ) .   Also  notice  that  the  first 
two  variables  entered  in  the  stepwise  discriminant  analysis 

(A  THK    and  AWET   )  were  subsequently  removed. 

The  importance  of  the  temperature  index  parameter  and  its 
time  differences  is  again  demonstrated  by  the  fact  that  they 
are  the  top  two  discriminators  in  these  equations.   The  mean 
values  for  ATI    were  about  1.4  C  for  foggers  and  0.6  C  for 
no-foggers.   Thus,  the  24-hour  change  in  the  temperature 
index  was  significantly  higher  for  fog  situations.   This 
result  is  likely  to  be  caused  by  increasing  temperatures  at 
the  top  of  the  inversion  which  is  associated  with  a  strength- 
ened inversion.   The  mean  values  for  A„TI^.  ^^  were  about 
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-5.2  C  and  -5.6  C  for  the  fog  and  no-fog  cases,  respectively. 
Although  the  BMD07M  program  found  some  statistical  signifi- 
cance in  the  group  means'  difference,  the  standard  deviations 
were  fairly  large  (about  2.4  for  both  groups)  and  very  few 
conclusions  can  be  made  about  this  parameter. 

The  950-mb  wind  direction  averages  were  about  214°  for 
fog  reports  and  186°  for  no-fog  cases.  Again  the  standard 
deviations  were  large  (nearly  115°  for  both  groups),  but  the 
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mean  wind  direction  prior  to  most  fog  occurrences  had  a  com- 
ponent from  the  west.  This  implies  that  the  fog  was  generally 
advected  landward  by  the  prevailing  southwesterly  flow. 

The  classification  functions  for  Quillayute's  wet  season, 
categories  1  and  2,  are  given  in  Table  XIX.   The  first  var- 
iable entered,  the  1600  PST  950-mb  wind  speed  (,SPD  ^),  was 
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the  best-discriminating  parameter  in  these  equations.   The 
mean  wind  speed  for  the  fog  cases  was  over  9.0  knots,  while 
the  no-fog  mean  was  5.7  knots.   The  cause  for  the  signifi- 
cantly higher  wind  speed  associated  with  fog  occurrences  is 
not  clear.   The  supposition  that  these  speeds  are  related  to 
pre-frontal  fog  appears  to  be  disproved  by  two  other  impor- 
tant discriminators. 

The  two-day-old  wind  direction  taken  at  0400  PST  (  DIR   ) 
was  the  second  best-discriminating  variable.   Surprisingly, 
the  mean  wind  directions  for  both  fog  and  no-fog  groups  were 
from  the  southeast  quadrant-- the  values  were  approximately 
155°  and  125°,  respectively.   The  third  important  variable 
was  the  sea-level  pressure  measured  at  1600  PST  (  SLP   )  . 
The  group  means  showed  a  significant  pressure  difference  of 
about  2.7  mb.   The  average  SLP  for  foggers  was  almost  1019.6 
mb ,  while  the  no-foggers  had  an  average  of  1016.9  mb .   This 
researcher  was  unable  to  find  an  adequate  explanation  for 
this  last  result  without  supplementary  information. 

Table  XX  shows  the  classification  functions  for  Quilla- 
yute  during  the  dry  season  for  categories  3  and  4.   The  F 
values  indicate  that  the  f irs t- entered  variable,  the  1600 
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PST  relative  humidity  (  RH   ) ,  is  by  far  the  most  important 
discriminator  in  these  equations.   The  mean  relative  humidi- 
ties for  the  fog  and  no-fog  cases  were  73.5%  and  60.4%,  res- 
pectively.  It  appears  logical  that  as  fog  approaches  the 
station  the  relative  humidity  will  increase.   The  13.1% 
spread  between  the  average  relative  humidities  provided  easy 
discrimination  between  the  two  groups. 

Three  other  variables  in  Table  XX  stand  out  as  signifi- 
cant in  the  classification  functions.   These  are  the  tempera- 
ture index  (.TI^^),  the  950-mb  wind  direction  (^DIR,^),  and 
lib  116 

the  moisture  index  ( ^ WET   ) ,  all  derived  from  the  most  recent 

1       ±  D 

1600  PST  observations.   The  TI  mean  values  were  -2.8  C  for 
foggers  and  about  -1.7  C  for  no-foggers.   The  fact  that  the 
fog  cases  had  a  larger  negative  TI  is  rather  surprising  since 
one  would  expect  that  the  closer  the  air  temperature  is  to 
the  SST,  the  greater  the  chance  of  fog.   The  negative  TI  in- 
dicates that  the  air  temperature  is  usually  lower  than  the 
SST. 

The  mean  wind  directions  for  the  fog  and  no-fog  groups 
were  both  from  the  southwest  (237°  for  foggers  and  231°  for 
no-foggers) ,  and  showed  such  a  minor  difference  that  one 
wonders  why  it  was  chosen  as  a  meaningful  discriminator.  The 
moisture  parameter  (  WET   )  seems  to  have  greater  theoretical 
meaning  in  this  program.   The  mean  moisture  values  were 
approximately  -8.3  C  for  fog  cases  and  -9.9  C  for  no-fog 
reports.   Recall  that  the  moisture  index  is  the  difference 
between  the  dew-point  and  sea-surface  temperatures.   In  the 
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Quillayute  dry  season  the  dew-point  temperature  is  nearly 
always  well  below  the  SST.   However,  during  fog  situations 
the  difference  is  smaller,  which  implies  a  greater  likeli- 
hood of  fog  formation. 

The  Quillayute  wet  season  classification  functions  for 
categories  3  and  4  are  presented  in  Table  XXI.   In  these 

equations  the  relative  humidity  parameter  ( ^ RH  ^)  again 

1   16 

stands  out  as  the  most  capable  discriminator.   As  before, 
the  difference  between  the  groups'  averages  is  about  13%, 
but  in  the  wet  season  the  relative  humidities  are  higher — 
almost  87.0%  for  foggers  and  74.2%  for  no-foggers. 

The  second-  and  third-best  discriminating  variables  are 
the  1600  PST  temperature  index  (  TI   )  and  the  0400  PST  wet- 
bulb  depression  ( ^ ^BD   ) .   The  TI  statistics  showed  fairly 
large  negative  mean  values  —  a  -7.7  C  for  fog  cases  and  a 
-8.1  C  for  no-fog  reports.   During  the  winter  season  one 
should  expect  larger  negative  values  for  this  parameter  in 
the  colder,  more  northerly  latitudes.   Nevertheless,  the 
slightly  smaller  TI  for  the  foggers  indicates  a  greater 
probability  of  a  fog  occurrence.   The  wet-bulb  depression 
means  values  showed  only  a  minor  spread  of  0.1  degree--the 
foggers  averaged  0.54  C  while  the  no-foggers  were  about 
0.64  C.   Although  the  difference  was  small,  the  discriminant 
analysis  routine  found  that  it  made  a  significant  contribu- 
tion in  these  classification  functions. 

Table  XXII  presents  the  Quillayute  fog  and  no-fog  dis- 
crimination percentages  by  season  and  category  pairs.   The 
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results  for  the  category  1  and  2  analyses  are  considerably 
better  than  for  3  and  4.   This  was  probably  due  to  the  much 
smaller  number  of  cases  in  the  1  and  2  category.   The  over- 
all weighted  discrimination  percentage  was  above  74.9  for 
the  fog  cases  and  nearly  75.5  for  the  no-fog  reports. 

The  cumulative  average  discrimination  percentages  for 
all  three  stations  analyzed  in  this  study  resulted  in  78.2% 
for  all  fog  cases  and  over  77.8%  for  all  no-fog  cases. 

E.    SOME  TEST  RESULTS 

In  order  to  test  the  discriminating  capability  of  some 
of  the  classification  equations  three  years  of  additional 
data  were  utilized.   The  Naval  Weather  Service  Detachment, 
Asheville,  North  Carolina  again  provided  the  airways  surface 
observations  and  radiosonde  soundings,  while  more  SST  infor- 
mation was  extracted  from  the  Fishing  Information  Bulletin 
Supplements.   The  test  data  covered  the  period  from  1  July 
1965  through  30  June  1968--i.e.,  the  three  years  preceding 
the  initial  analysis  period. 

Primarily  due  to  time  considerations,  only  the  San  Diego 
classification  equations  were  tested.   The  data  were  prepared 
by  season  and  categories  in  the  manner  previously  described 
(see  Section  IV,  subsection  A).   The  number  of  fog  and  no-fog 
cases  listed  by  category  for  the  test  data  dry  and  wet  sea- 
sons are  shown  in  Tables  XXIII  and  XXIV,  respectively.   A 
testing  computer  program  was  developed  which  would  form  the 
fog  and  no-fog  probabilities  (see  Eqs.  (1)  and  (2))  and 
classify  each  case  in  the  same  manner  as  the  BMD07M  routine. 
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Table  XXV  shows  the  percentage  of  correct  discriminations 
which  resulted  from  testing  the  San  Deigo  classification  equa- 
tions on  the  independent  data  set.   Although  the  discrimina- 
tion percentages  were  considerably  reduced  from  the  dependent 
analysis  sample,  the  outcome  was  generally  acceptable  except 
for  the  dry  season  fog  prediction  capability.   The  reason  for 
this  failure  is  not  known.   The  overall  weighted  fog  discrim- 
ination percentage  for  the  San  Diego  test  data  was  above 
58.4,  while  the  no-fog  percentage  was  nearly  74.3. 

By  comparing  the  numbers  of  fog  and  no-fog  cases  in 
Tables  III  and  XXIV,  one  will  note  that  the  relative  propor- 
tions of  cases  are  similar.   Thus,  it  seems  unlikely  that  un- 
usual group  sizes  could  be  a  factor  in  the  dry  season's  poor 
discrimination  capability.   Perhaps  better  results  could  be 
obtained  if  the  BMD07M  program  prior  probabilities  were  made 
proportional  to  the  number  of  cases  in  each  group.   This 
would  increase  the  likelihood  of  being  assigned  to  the  larger 
group.   Another  alternative  would  be  to  make  a  Bayesian  ad- 
justment to  the  probabilities  of  group  membership,  especially 
when  the  groups  are  of  grossly  different  sizes  as  was  found 
in  both  San  Diego  and  Oakland. 

From  Table  XXIII  one  can  find  that  during  San  Diego's 
dry  season  only  11%  of  the  total  number  of  cases  (for  the 
test  data  period)  were  fog  reports.   During  the  wet  season 
(see  Table  XXIV)  about  twice  as  many  (22.5%)  fog  days  oc- 
curred.  From  these  statistics  a  forecaster  might  logically 
conclude  that  he  would  be  correct  almost  90%  of  the  time  in 
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the  dry  season  by  never  predicting  the  occurrence  of  fog! 
Indeed,  persistence  and  climatology  forecast  methods  may 
be  difficult  to  surpass,  but  other  prediction  tools  are 
essential . 
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VI .    PROCEDURE  FOR  APPLICATION 
OF  THE 
CLASSIFICATION  FUNCTIONS 


One  of  the  primary  objectives  of  this  study  was  to  de- 
velop classification  equations  which  could  be  utilized  to 
forecast  the  daily  probability  of  occurrence  of  marine  fog 
at  the  three  U.  S.  west  coast  stations  analyzed.   The  clas- 
sification functions  presented  herein  provide  a  means  for 
forecasting  the  probability  of  fog/no  fog  for  a  24-hour 
period  beginning  at  1601  PST.   In  order  to  calculate  these 
probabilities  the  forecaster  must  have: 

1.  The  station's  3-hourly  surface  observations  and  the 
upper-air  soundings  for  the  last  three  days. 

2.  A  reliable  bi-weekly  mean  coastal  SST  near  the  sta- 
tion.   (This  is  required  to  form  the  TI  and  WET  parameters.) 

3.  A  desk  calculator  or  a  small  programmable  computer. 
The  following  step-by-step  procedure  is  recommended  for 

the  proper  application  of  the  classification  functions. 

1.  Note  whether  fog  has  occurred  on  any  surface  report 
in  the  last  three  days,  starting  with  the  most  recent  1600 
PST  observation.   If  there  has  been  such  a  fog  occurrence, 
use  the  table  for  the  current  season  (wet  or  dry) ,  cate- 
gories 3  and  4.   If  there  is  no^  history  of  fog,  use  the 
table  labeled  categories  1  and  2,  for  the  current  season. 

2.  Using  the  surface  observations,  the  upper-air  sound- 
ings, and  the  SST  data,  find  the  numeric  value  of  each  of 
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the  variables  listed  in  the  table  selected  in  step  1.  (See 
the  definitions  of  the  variables  in  Section  IV,  subsection 
B.) 

3.  With  the  aid  of  a  desk  calculator  or  a  computer  mul- 
tiply the  raw  values  of  the  variables  by  their  associated 
fog  coefficients  and  add  them  together  along  with  the  con- 
stant.  Do  the  same  with  the  no-fog  coefficients.   These  two 

sums  are  the  discriminant  scores,  D^  and  D  ^  (see  Eq.  (1)) . 

f        nf 

4.  To  avoid  problems  in  exponentiating  discriminant 

4 
scores  which  might  be  outside  the  allowable  range,   one  may 

do  the  following: 

a.  Select  the  minimum  discriminant  score  between 

D^  and  D  ^. 
f        nf 

b.  Subtract  the  minimum  discriminant  score  from  D  . 
Let  the  result  be  called  "DIFF". 

(1)  If  DIFF  is  greater  than  or  equal  to 
-174.673,  exponentiate  DIFF. 

(2)  Otherwise,  set  DIFF  equal  to  zero  and  then 
exponentiate . 

Call  the  above  result  "PF". 

c.  Repeat  step  b  using  the  no-fog  discriminant  score, 

D  ^,  in  place  of  D^ .   Call  this  result  "PNF". 
nf  f 

5.  Form  the  posterior  probability  of  a  fog  occurrence 
by  dividing  PF  by  the  sum  of  PF  and  PNF  (see  Eq .  (2)).   Then 


The  constant  174.673  is  the  largest  power  of  e  that 
the  IBM  360  computer  can  represent.   If  using  another  com- 
puter this  constant  will  vary. 
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the  probability  of  no  fog  occurring  is  unity  minus  the  , pro- 
bability of  fog. 

Appendix  A  presents  a  generalized  FORTRAN  program  of 
the  above  procedure. 
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VII .    CONCLUSIONS  AND  FINAL  REMARKS 

The  primary  goal  of  developing  linear  classification 
equations  which  may  be  used  for  forecasting  the  probability 
of  marine  fog  at  San  Diego,  Oakland,  or  Quillayute,  has  been 
accomplished.   The  ultimate  capability  of  these  equations 
has  not  been  fully  tested.   Nevertheless,  it  is  hoped  that 
these  discriminant  functions  may  serve  some  useful  purpose — 
if  not  for  operational  employment,  then  perhaps  as  a  basis 
for  further  study. 

One  of  the  important  by-products  of  the  classification 
equations  v/as  insight  as  to  which  variables  are  the  best 
fog/no-fog  discriminators.   The  following  is  a  brief  summary 
of  the  most  capable  discriminating  variables,  listed  in 
order  of  decreasing  importance  by  location,  season,  and 
category  pair: 
SAN  DIEGO 

DRY  (1  +  2):    ^TI^^     ^THK^^     ^2^^^6-04     ^^^%4 
WET  (1+2):    ^WBD^^   3HT^^      ^TI^^  ATB^^ 

DRY   (3  +  4):    ^TI^^     3°^^04     ^^WB^^^         ^^BD^^ 
WET  (3+4):    ^WBD^^   ^HB^^      ^FS^^  ATI^^ 

OAKLAND 

DRY   (1+2):    A^THK^^  A^RH^^     A^TI^^         ^2^^04-16 
WET  (1  +  2):    ^TT^g    ^RH^^      ^^%4-16      3°^^6 
DRY   (3  +  4):    ^FS^^     ^^PD^^     ^TWB^^         ^^^16-04 
WET  (3  +  4):    ^TT^^    ASPD^^    ^HT^^  ^^^^04-16 
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QUILLAYUTE 

DRY   (1+2) : 
WET  (1+2) : 

DRY   (3  +  4)  : 
WET  (3+4) : 


^^^04   ^2^^04-16  1^^^04 


l^^^e   2^^%4 


l^"l6 


TI 
1   16 


i^^6   I'^^ie 


2^^^6 


1°^^6 


1     04 


'^""16 


Several  conclusions  can  be  made  about  the  discriminating 
parameters  shown  in  this  analysis  summary: 

1.  The  temperature  index  (TI)  and  several  of  its  time 
differences  are  powerful  discriminators,  especially  in  San 
Diego  and  Quillayute„   Since  this  index  uses  the  SST  as  an 
"anchor  point",  reliable  coastal  SST  data  are  very  important. 
This  may  provide  a  clue  as  to  why  the  TI  is  not  particularly 
significant  at  Oakland.   The  SST  data  for  this  station  was 
extracted  from  SST  charts  at  a  point  just  west  of  the  San 
Francisco  Bay  entrance,  supposedly  where  the  marine  fog 
might  form  and  be  advected  shoreward.   In  retrospect,  per- 
haps a  bay  water  temperature  would  have  been  a  better  refer- 
ence for  the  TI  at  this  station. 

2.  The  most  recent  relative  humidity  measurement  ( ^ RH ^ ^ ) 

115 

and  its  24-hour  change  (A  RH^.)  appear  to  be  excellent  dis- 

1   Id 

criminating  variables  in  Quillayute  for  categories  3  and  4, 
and  in  Oakland  for  categories  1  and  2. 

3.  In  Oakland's  wet  season,  the  1600  PST  temperature 
at  the  inverstion  top  (  TT  )  is  a  useful  fog/no-fog  pre- 
dictor . 

4.  The  1600  PST  wet-bulb  depression  is  an  important  dis- 
criminator in  San  Diego  during  the  wet  season. 
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5.   Two  discriminating  variables  which  were  hoped  to 
be  of  importance  in  the  final  classification  equations  were 
the  strength  of  the  inversion  (SI)  and  the  moisture  index 
(V7ET)  .   One  will  notice,  however,  that  these  parameters 
generally  made  only  a  relatively  minor  contribution  to  the 
discriminating  capability  of  most  equations.   Several  times 
the  moisture  index  was  entered  into  the  regression  equation 
first,  indicating  a  large  initial  discrimination  power,  only 
to  loose  its  importance  as  other  variables  were  entered. 
Nevertheless,  this  researcher  feels  that  the  moisture  index 
is  still  a  worthwhile  parameter  and  should  be  investigated 
in  further  research. 

The  discriminant  analysis  approach  is  certainly  not  a 
new  research  technique.   In  the  past,  multiple  linear  regres- 
sion has  been  employed  in  numerous  studies  which  require 
analysis  of  multivariate  observations  whose  predictand  was 
non-numerical.   In  the  last  decade,  through  the  advance  of 
computer  technology,  tremendous  progress  has  been  made  in 
the  rapid  analysis  of  numerical  and  statistical  information. 
Without  the  use  of  a  high-speed  computer  and  a  sophisticated 
discriminant  analysis  program  the  completion  of  this  re- 
search by  a  single  individual  within  a  reasonable  time  frame 
would  be  impossible.   However,  through  the  application  of 
these  powerful  computational  resources  meaningful  relation- 
ships can  be  determined  and  an  enhanced  f og- f or ecas t i ng 
capability  may  be  evolved. 
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VIII.    RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  analysis  results  revealed  several  interesting  and 
unexpected  relations  between  various  meteorological  para- 
meters and  the  occurrence  or  non-occurrence  of  marine  fog 
at  selected  locations  on  the  west  coast  of  the  United  States 
In  order  to  evaluate  the  full  significance  of  these  para- 
meters and  attempt  to  improve  the  classification  functions, 
the  following  recommendations  are  offered  for  future  study. 

1.  A  larger  data  base  is  needed  in  order  to  obtain  a 
greater  number  of  fog  and  no-fog  cases  for  discriminant 
analysis.   A  distinct  paucity  of  fog  cases  was  noted  in 
both  San  Diego  and  Oakland  dry  seasons.   In  such  instances 
the  BMD07M  program  may  generate  classification  equations 
which  have  used  a  few  anomalous  occurrences  (which,  however, 
may  be  a  significant  portion  of  a  small  data  set)  for  dis- 
crimination purposes.   Such  equations  might  not  indicate 
the  discriminating  variables  which  are  truly  important  to 
fog  forecasting.   Consequently,  a  ten-year  data  period  is 
considered  to  be  the  minimum  for  best  results  using  the  dis- 
criminant analysis  technique. 

2.  Make  adjustments  in  the  prior  probabilities  such 
that  these  probabilities  are  proportional  to  the  number  of 
cases  in  each  group.   Further  experimentation  with  variable 
prior  probabilities  is  necessary  to  test  their  effect  on 
the  selection  of  the  fog/no-fog  discrimination  parameters 
and  the  final  classification  equations. 
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3.  Expand  the  list  of  discriminating  variables  to  in- 
clude visibility,  low  cloud  information  and  persistence 
parameters.   Other  wind  levels  should  also  be  studied. 

4.  Refinement  or  re-definition  of  "fog"  to  include 
only  those  cases  where  visibilities  are  reduced  to  less 
than  three  miles  (IFR  conditions)  may  result  in  forecasts 
of  greater  operational  importance.   Some  measure  of  rela- 
tive fog  intensity,  such  as  "heavy"  versus  "light"  fog  may 
also  enhance  the  forecast.   However,  a  very  large  data 
base  will  be  required  for  such  refinements. 

5.  The  use  of  better  sea-surface  temperature  informa- 
tion is  strongly  recommended.   Bi-weekly  average  SST  data 
extracted  from  smoothed  isotherm  analyses  are  likely  to  be 
inaccurate.   The  SST  data  are  important  because  a  one-  or 
two-degree  temperature  difference  between  the  marine  air 
and  the  underlying  sea  may  be  crucial  in  determining 
whether  fog  is  expected  to  form,  or  in  measuring  its  inten- 
sity.  Reliable,  dai ly  SST  information  might  make  the  tem- 
perature and  moisture  indices  (TI  and  WET)  better  fog/ 
no-fog  discriminators. 

6.  With  the  use  of  a  larger  data  set,  one  should  try 
analysis  on  four  seasons  of  the  year  instead  of  just  two. 
Such  a  seasonal  refinement  may  provide  an  enhanced  fog  pre- 
diction capability  through  the  use  of  more  representative 
classification  functions. 

7.  Develop  a  method  of  stratifying  the  24-hour  fore- 
cast into  two  12-hour  periods  (or  smaller  intervals)  in  an 
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effort  to  indicate  more  precisely  the  expected  time  of  fog 
occurrence . 

8,   Since  Oakland  is  somewhat  sheltered  from  the  marine- 
air  influence  by  the  coastal  hills  of  San  Francisco,  data 
from  a  station  closer  to  the  coast  should  be  used.   The 
surface  observations  from  such  a  station  could  be  coupled 
with  the  Oakland  upper-air  sounding  data. 

Application  of  these  refined  analysis  techniques  to  an 
enlarged  data  base  would,  hopefully,  result  in  a  significant 
improvement  in  marine-fog  forecasting  at  the  west  coast 
stations  of  interest. 
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CATEGORY 

HISTORY  (3  DAYS) 

CURRENT  24-nR  PERIOD 

1 
2 

NO  FOG 
NO  FOG 

FOG 
NO  FOG 

3 

4 

FOG 
FOG 

FOG 
NO  FOG 

TABLE  I.    FOG  AND  NO-FOG  CATEGORIES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 

50  FOGGERS 

2 

531  NO-FOGGERS 

3 

60  FOGGERS 

4 

180  NO-FOGGERS 

TABLE  II.    SAN  DIEGO  DRY  SEASON  FOG  AND  NO-FOG  CASES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

69  FOGGERS 
407  NO-FOGGERS 

3 
4 

108  FOGGERS 
220  NO-FOGGERS 

TABLE  III.    SAN  DIEGO  WET  SEASON  FOG  AND  NO-FOG  CASES 
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VARIABLE 


Atvjb 


16 


l"04 


AWED 

Adpd 

ASLP 


04 
04 


16 


A^TWB^ , 
2  04 

2«'^04 

2^^04 

A, TDP, ,     ^, 
1  16-04 

A    TI 
1        16 

TB 
116 

1^^^04 

1^^^6 

-■  TWB^  , 
1  04 

Constant 


Coefficients  of  the 


FOG  FUNCTION 


-0.69860 

-0. 73689 

0. 01075 

0. 00125 

-1. 23321 

0. 34705 

-0. 36144 

0.11503 

0. 00153 

0. 01062 

1. 19864 

-0.41794 

-0.24155 

-0. 00005 

0.14416 

2.35226 

37.83261 


NO-FOG  FUNCTION 


-0. 94497 

-1.07313 

0.00646 

0. 00313 

0.49745 

-0.27131 

-0.17900 

-0. 09174 

0.00081 

0. 11482 

1.00863 

-0. 54457 

-0. 16593 

0.00234 

0. 13815 

2  .  24390 

-35.13945 


FINAL  F  VALUE 


2.7169 
23.6201 
9.  2365 
10. 2577 
8.9386 
5.1529 
8. 5494 
1.9007 
3  .8738 
3.7899 
2.9281 
1.8798 
4 .2166 
2.8459 
2. 2897 
1.8560 


TABLE  IV.    SAN  DIEGO  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFICA- 
TION FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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VARIABLE 


2  =  ^04 

A  TT 
2   04-16 

1     04 
3"%4 


Ate 


04 


3^^04 

3=^°16 

1  04 

A^HB, ,     ^  , 
2        16-04 

2  04 
3^16 


ASI 


16 


2^04 

2  =  ^04 
Constant 


FOG 


Coefficients 
FUNCTION   f  NO- 


of  the: 
•FOG  FUNCTION 


FINAL  F  VALUE 


(Removed  at 
(Removed  at 

1. 22214 

0.00375 
-lo64342 

(Removed  at 

2.40756 
-1. 13487 

0. 02395 

0. 00008 
-0. 07748 

0.60771 
-0.47578 

(Removed  at 

2. 17704 

1.17041 

0.01875 

0. 00054 

(Removed  at 

0.19322 
-2.96269 

1 .05731 
-3. 10002 

0.06976 
-1.85104 

1.48530 
■40. 16452 


step  16,  &  later  re-entered) 
step  20,  &  later  re-entered) 


0.99060 

0.00260 

-1.27683 

step  22) 
2. 55562 

-1  .  29938 
0.02075 
0.00134 

-0.19115 
0.49668 

-0.64106 

step  25 ) 
1.65256 
1.27535 
0.01652 
0.00094 

step  31) 
0. 16090 

-1. 90545 
0.98900 

-2. 83606 
0.19479 

-1.68469 
1.39985 

40.49783 


2  .9369 
4.6784 
7 . 7045 

3 .4998 
4.0756 
3.4434 
10. 5531 
7.8624 
4 .1412 
4 .0329 

9.7573 
1.9934 
2  .4990 
3.  2675 

1. 5050 
12. 1110 
2.0290 
3.0366 
1. 5140 
3.1752 
1.6265 


TABLE  V.    SAN  DIEGO  V?ET  SEASON  FOG  AND  NO-FOG  CLASSIFICA- 
TION FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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VARIABLE 


1^^6 

1   04 

1  04 

TT 

2  04 

AWBD^ . 
04 


Coefficients  of  the 


Arh 


04 


^^^^6-04 
^2^^^04-16 

^"^6-04 

2^^^6 
1     04 


Atwb 


16 


A  TDB 

2  15 

3"l6 
l''^%4 
'^2''%4 

3  =  ^04 


ATI 

Ahb 


04 


16 
Constant 


FOG  FUNCTION 


127. 58673 
0.42659 


NO-FOG  FUNCTION 


FINAL  F  VALUE 


126.82245 
0.43420 


(Removed  at  step  21) 


-23.87572 
-7. 25344 
26. 61739 
32.45081 
10. 95739 

-38.42950 
24.53719 


-24.84428 
-8. 02159 
26.73076 
30.04062 
10. 76887 

-38. 25275 
24. 32675 


(Removed  at  step  19) 

-0.01331    I      -0.01411 

(Removed  at  step  25) 


75.30409 

0.41899 

221. 02830 

-1»12342 

-97 . 16171 

39. 62401 

68. 23961 

59.41071 

17.00053 

-17. 24670 

-5. 97035 

-24. 14389 

-51.37943 

0.01494 

112260. 12500 


75. 55360 

0.42067 

221 . 20645 

-1. 12795 
-97.94383 
40. 02596 
68.77629 
59.89287 
17.34978 
-17. 02985 
-6.00543 
-23 . 99202 
-51. 23575 
0.01441 
-112446.18750 


3 .4129 
19.8636 

22 . 1792 
7. 1250 
6.  2731 

10. 2524 
4  .  3374 
2. 0741 
3 .0246 

1.7733 

4.1700 
3 .4352 

3  .4470 

1.4305 
12. 1096 

6. 3403 
10. 1973 

7. 0221 

4  .  3841 
1.7932 
1. 2460 
2.  3931 
1.3127 
1.1024 


TABLE  VI.    SAN  DIEGO  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFICA- 
TION FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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Coefficients  of  the: 

VARIABLE 

FOG  FUNCTION    NO-FOG  FUNCTION 

FINAL  F  VALUE 

l""^16 

(Removed  at  step  23) 

1^=16 

0.27487 

-0.50531 

10.3541 

i^^e 

-0.04832 

-0.04735 

11.0117 

'""is 

-25.14560 

-25. 54906 

8.6549 

1^"'^04 

0.15735 

0.15587 

3  .6692 

A„TWB^  . 
2     16 

(Removeca  at  step  22) 

ASPD, ^ 

1  D 

16.75336 

16.83115 

3.5183 

^2"%4 

0.01693 

0.01492 

3  .  8820 

^  =  ^16 

(Removec3  at  step  19) 

1°^''04 

-0.14458 

-0.14257 

1.9764 

Atwb, ^ 

1 6 

-1.23822 

-0.85878 

4 .6643 

'^''"16 

1.19872 

1. 13272 

7.7717 

A^TDB^^ 

20.99492 

20.64488 

6.1341 

^2''«04 

-2.62901 

-2.64525 

2.2931 

i^'^^e 

112.41910 

112.31770 

3.9060 

^^^04 

-5.89893 

-5.82536 

3 . 2766 

^^2^04 

0.90936 

0.79580 

3 .9585 

"2^^16 

-6.86152 

-6.61983 

2.9379 

2™''l6 

53. 50967 

53.32607 

4 .8779 

I'^'^e 

20. 14551 

20.63298 

14.  2448 

I'^^ie 

-7.69770 

-7. 90987 

4.3569 

3^^6 

-19.78639 

-19.68510 

1.0576 

'^2^^16-04 

-1.14815 

-1.20561 

2.7393 

*2«^6-04 

-0.01367 

-0.01401 

1. 5600 

2^^^6 

10.54596 

10.48309 

1.3276 

Constant  -I 

>7540. 84766 

-57436.34766 

TABLE  VII. 


SAN  DIEGO  WET  SEASON  FOG  AND  NO-FOG  CLASSIFICA- 
TION FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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SEASON  &  CATEGORY 

FOG  PERCENTAGE 

NO-FOG  PERCENTAGE 

DRY   (1+2) 
WET  (1+2) 
DRY   (3+4) 
WET  (3+4) 

78.0 
75.4 
85.0 
87.  0 

79.3 
81.6 
86.1 
81.4 

WEIGHTED  AVERAGE 

82.2 

81.3 

TABLE  VIII. 


SAN  DIEGO  FOG  AND  NO-FOG  DISCRIMINATION 
PERCENTAGES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

25  FOGGERS 
692  NO-FOGGERS 

3 

4 

7  FOGGERS 
78  NO-FOGGERS 

TABLE  IX.    OAKLAND  DRY  SEASON  FOG  AND  NO-FOG  CASES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

66  FOGGERS 
345  NO-FOGGERS 

3 

4 

123  FOGGERS 
246  NO-FOGGERS 

TABLE  X.    OAKLAND  WET  SEASON  FOG  AND  NO-FOG  CASES 
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Coefficients  of  the: 

VARIABLE 

FOG  FUNCTION 

NO-FOG  FUNCTION 

FINAL  F  VALUE 

^l^"l6 

0.12679 

0.06999 

6.8501 

A^THK^^ 

-0.00299 

-0.00066 

9.4800 

2"^^6 

1  .  05340 

1.28505 

4.5235 

^2^^^04 

-0.13446 

0.01399 

3. 5524 

^TDB^^ 

-0.22813 

0.06931 

3.9888 

^^^04 

(Removed  at  step  15) 

^2^^04-16 

0.07802 

0.22367 

4.9993 

^^^o4 

(Removec 

i  at  step  16) 

1^^04 

0. 59826 

0.52878 

2.1779 

^2^^^6 

0.00227 

-0.00105 

1.1546 

2^^^04 

0.04460 

0.04094 

1.8274 

3"^6 

0.00316 

0.00258 

1.6784 

3^^%4 

0.63973 

0.76082 

1.9224 

^2^^04 

0.08523 

-0.07332 

5.2857 

Constant 

-13 .89205 

-13.44692 

TABLE  XI.   OAKLAND  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFICATION 
FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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Coefficients  of  the: 

VARIABLE 

FOG  FUNCTION 

NO-FOG  FUNCTION 

FINAL  F  VALUE 

l^"l6 

3.34849 

3.27486 

10. 2524 

i^^ie 

-3  .  08567 

-3.17853 

16.3925 

1^^^04 

42. 07230 

42.00005 

5.6511 

^^^04 

-9.95635 

-9.84301 

5.4529 

3^'^16 

-Ool7242 

-0. 16732 

6.7092 

^^%4-16 

0.02788 

-0o04618 

7.7952 

A^WET^^ 

-4 .01388 

-4.08238 

2.8188 

^2^^^6 

-0. 11382 

-O0IIO6O 

4.4790 

^^^%4 

(Removed 

[  at  step  17) 

^2^^04-16 

8.11221 

8. 20044 

2.  5340 

^"^6 

-0.00199 

-0. 00223 

2.8280 

A^THK^^ 

-0. 09082 

-0.09181 

2.5174 

3^^^04 

0. 13206 

0.13412 

1.9915 

l"^04 

0.00785 

0.00697 

4.7836 

^"%4 

0.02822 

0.02869 

3.0078 

^^^^6-04 

7.96099 

7.82930 

2.5799 

3^"^04 

-0.03665 

-0.03714 

1.  2301 

1^^^04 

9o87953 

9o94415 

1.1570 

Constant 

-21517.31641 

-21437.86719 

TABLE  XII. 


OAKLAND  VJET  SEASON  FOG  AND  NO-FOG  CLASSIFICATION 
FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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Coefficients  of  the: 

VARIABLE 

FOG  FUNCTION 

NO-FOG  FUNCTION 

FINAL  F  VALUE 

^■'P%4 

4.41206 

8.59220 

5.7808 

AWBDp^ 

-9.61003 

-15.05353 

2.6278 

1«^16 

0. 00639 

0.00851 

2.2578 

'^2'^^^6 

(Removed  at  step  12) 

^"l6-04 

-0.19130 

0.27583 

6.8050 

2''%4 

-1.24039 

-0.81018 

2.3739 

2^^16 

1„40700 

0.73658 

6.6531 

^^^%4-16 

(Removed  at  step  18) 

2''^6 

-0. 96936 

3. 33698 

12.4420 

2™^6 

7.61710 

5.66939 

8.4864 

'^2^^^4 

-0.30712 

0.30540 

2.9913 

*2^^%4-16 

(Removec 

i  at  step  21) 

1™^6 

2.40190 

3  .20982 

3.0565 

^2^^04 

-0.42721 

-0.65045 

1.7665 

2^^''l6 

lo35426 

2.62428 

9.8126 

I'^^e 

-1.02120 

-0.61391 

5.4850 

'^I^IS 

0.44584 

-0cl6509 

2.7886 

'^2  =  ^%4 

-0.04702 

-0.60388 

3.4339 

Constant 

-76.13808 

-73.33984 

TABLE  XIII 


OAKLAND  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFICATION 
FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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VARIABLE 


Coefficients  of  the: 
FOG  FUNCTION     NO-FOG  FUNCTION 


FINAL  F  VALUE 


l^"l6 

ASPD, ^ 
15 

A, SLP^  ,     ,  ^ 
1     04-16 

3='=04 

A^TT^ , 
2       04 

''2^^^4 


Adir 


04 


^2''^''04-16 
'^2°'''=16 
'^2''«16 

1™^4 

2^"04 
Constant 


22.  77654 

-11 .35367 

2.84706 

-5.79911 

0. 23442 

31. 59393 

-0.01989 

64. 52661 

-2.40382 

(Removed 

0.02141 

7.10984 

-7. 96592 

-4.41379 

-32.83693 

-103.94606 

-57o20595 

0. 38270 

0. 21845 

-0. 53928 

-1.36683 

21.05357 

-19.81030 

69. 76450 

38127.37891 


22. 63602 

-11.43764 

2. 60788 

-5. 69068 

0. 23942 

31. 41635 

-0.01919 

64.19128 

-2.48583 

at  step  23 ) 

0. 02209 

6.99197 

-8.  12432 

-4.44339 

-32. 94241 

-104.36403 

-57.47128 

0.38461 

0. 21671 

-0.87731 

-1.42692 

21 .17676 

-19.91208 

69.88263 

37993. 56641 


7 . 2474 
11.3471 
5.6923 
10.4081 
6.5569 
8.6105 
9.8690 
2.6085 
8.0632 

3 .4260 

2.4149 
1.0617 
1.0977 
2.7667 
2.9697 
2 . 5596 
1.6718 
2 . 0633 
2.7606 
1.4863 
4.4619 
1.7678 
1.2062 


TABLE  XIV. 


OAKLAND  WET  SEASON  FOG  AND  NO-FOG  CLASSIFICATION 
FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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SEASON  &  CATEGORY 

FOG  PERCENTAGE 

NO-FOG  PERCENTAGE 

DRY   (1+2) 
WET  (1+2) 
DRY   (3+4) 
WET  (3+4) 

76.0 
78.8 
85.  7 
76.4 

74.6 
73  .0 
98.7 
81.3 

WEIGHTED  AVERAGE 

77.4 

76.8 

TABLE  XV.    OAKLAND  FOG  AND  NO-FOG  DISCRIMINATION  PERCENTAGES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

47  FOGGERS 

67  NO-FOGGERS 

3 

4 

290  FOGGERS 
270  NO-FOGGERS 

TABLE  XVI.    QUILLAYUTE  DRY  SEASON  FOG  AND  NO-FOG  CASES 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

39  FOGGERS 

91  NO-FOGGERS 

3 

4 

327  FOGGERS 
225  NO-FOGGERS 

TABLE  XVII.    QUILLAYUTE  WET  SEASON  FOG  AND  NO-FOG  CASES 
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VARIABLE 


Coef f i  cien ts 
FOG  FUNCTION     NO 


of  the: 
-FOG  FUNCTION 


FINAL  F  VALUE 


^2^'"'^04 


AWET 


16 


'^2^^04 
3"%4 
3  =  ^6 

3"^6 

A,TI„^ 

1  04 

'^2'*"l6 
^  =  ^■^16 
*2^^'^04 
2"^%4 
^^%4 

^2^^04-16 

SLP 

2  16 

A  TWB 
116-04 

ARH 

^^^%4 
ATDP^g 

^"l6-04 

^^"^04 
^2^^16-04 

l"%4 
^2^^04-16 

1°^%4 

2»^6 
Constant 


(Removed  at 
(Removed  at 
-12.59450     I 

(Removed  at 

0.02701 

23 .46292 

0. 05579 

-54. 50514 

-2.06415 

0.02754 

23.25471 

9.75030 

(Removed  at 

(Removed  at 

-58.78047 

85.17726 

-21.78006 

-9.34072 

1.33832 

27 . 53899 

-50.82481 

-0. 20962 

0.76721 

0.13531 

-0.01931 

-20c46352 

-Oc  09945 

43064.78516      -43 


step  28) 
step  18) 

-12.84844 

step  23) 

0. 02604 

23 . 03337 

0. 05850 

-55. 63574 
-2. 08957 
-0. 03313 
23 . 53885 
10. 26355 
step' 25) 
step  29) 

-60. 06639 
85.38756 

-22. 06129 

-9. 59527 

1.35184 

28. 27422 

-51o49023 

-0. 21217 

0.85447 

0.13775 

-0. 01804 

-20.74287 
-0.09997 

273.94141 


3. 5198 

3 .1081 

2.9156 

5. 3234 

29.4205 

21.8740 

4.8483 

5.7243 

10.0187 


23  .  5783 

6. 3899 

1.9874 

16.8019 

10.7870 

11 .6391 

7  .8411 

4 .7505 

2.8860 

8. 0030 

1.7703 

1. 3205 

1.0058 


TABLE  XVIII. 


QUILLAYUTE  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFI- 
CATION FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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VARIABLE 


2"^16 
2°^''04 

1   04 

-,  THK^  , 
1     04 

-,DIR^  , 
3     04 

^^2^^6-04 

1^^%4 

'^2^04-16 

^2=f°16-04 

^■^^^6-04 

^^2=^16 
Constant 


Coefficients  of  the: 


FOG  FUNCTION 


NO-FOG  FUNCTION 


13.96301  13.56759 

18.23346  18.10907 

1.03278  0.90201 

-0.07289  -0.07226 

(Removed  at  step  21) 

(Removed  at  step  16) 


0. 18567 


0. 16958 


(Removed  at  step  14) 


-21. 14665 


-20. 92448 


(Removed  at  step  19) 


-0. 20689 

0. 10502 

-0.10421 

0o38224 

-4.36938 

13 . 52794 

2.10825 

-5. 01223 

-4  .28529 

-9291.10938 


-0. 19983 

0. 11072 

-0.10166 

0. 25265 

-4.17476 

13  .18516 

2.01247 

-5. 12174 

-4.10700 

-9163. 24219 


FINAL  F  VALUE 


27.1404 

14.1167 

2. 4866 

2 . 3006 


14.3619 

3.6018 

4.  2505 
4. 7000 
5.5468 
6.  2855 
4.  9473 
5.4682 
1.9396 
1 . 3300 
1.1963 


TABLE  XIX.    QUILLAYUTE  WET  SEASON  FOG  AND  NO-FOG  CLASSIFI- 
CATION FUNCTIONS  FOR  CATEGORIES  1  AND  2 
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VARIABLE 


l^"l6 


ATI 


04 


^2^^04 

A  FS 
2   16 

A  TDP 
2     04-16 

A  THK^^ 
1     04 


Adpd 


04 


1^^^6 
^2^^^04-16 
\^^^6-04 


Ahb 


16 


3™^6 

^™%4-16 
1^^°16 


Ahb 


04 


,HB^  ^ 
1   04 

Cons  tant 


Coefficients  of  th( 


FOG  FUNCTION 


10.36016 
3.93438 

26.64711 
-21.  50443 

-0. 62858 
0.10819 

-0.81858 
1.16413 

-6.98908 

-0. 00598 
0.47370 

49. 51489 
1.84741 
0. 62066 

-0. 01738 
0. 02608 
7.  12852 

-0.75415 

-0.49160 

14.06279 

-7.78271 
0. 01489 

-0.01064 
25618.90625 


NO-FOG  FUNCTION 


10.21278 
3.88034 

26. 28835 
-21.  19310 

-0.63444 
0.11094 

-0.87311 
1.00654 

-6.84136 

-0.00620 
0. 63786 

49. 56006 
1.80900 
0.62268 

-0. 01752 
0.02658 
7. 20581 

-0.79099 

-0.63447 

14.09487 

-7.70232 
0.01461 

-0. 01032 
25653 .00781 


FINAL  F  VALUE 


34.0688 
2.8793 

12.  3250 
8.5402 
9.6180 
5.5265 
8  .3780 
4.2382 
7.3719 
2.6089 
3  .9646 
4.1959 
2.2311 
1„5758 
1.9450 
1.4060 
3.  2643 
1.9279 
3 . 5914 
1.1353 
1.9297 
2. 7728 
1 .8722 


TABLE  XX.    QUILLAYUTE  DRY  SEASON  FOG  AND  NO-FOG  CLASSIFI- 
CATION FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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Coefficients  of  the: 

VARIABLE 

FOG  FUNCTION 

NO-FOG  FUNCTION 

FINAL  F  VALUE 

i-'Nie 

4.04704 

3.95495 

75.9569 

I'^^e 

-6.03803 

-6. 21370 

19.3615 

3^^6 

2.93756 

2.72885 

7.0713 

1^^6 

(Removed  at  step  15) 

^■^"04 

-0.12554 

-0. 12898 

7.6744 

^"%4-16 

0.01320 

0.01337 

2.0153 

A, SPD,  ^  ^  ^ 
1     16-04 

-1.65702 

-1.71154 

7.9495 

l"^%4 

30.97249 

30.44763 

10.2992 

1^^4 

-1.81256 

-2.01150 

5.9802 

3^"^04 

-0. 02820 

-0.02772 

4.1084 

^°^^16 

0.  20934 

0. 20692 

5.8071 

1^^^04 

13  .66911 

13.64324 

4.7360 

3^^16 

1.19262 

1.13235 

5.8522 

3^^^6 

0.04358 

0.12552 

4.2738 

^1°'^04 

0.  02501 

0.  02693 

1.9975 

^"%4 

-0. 00430 

-0. 00412 

1.9398 

'^i^^^ie 

-2.32839 

-2.31851 

1.2658 

Constant 

-7134.03516 

-7101. 53516 

TABLE  XXI.    OUILLAYUTE  WET  SEASON  FOG  AND  NO-FOG  CLASSIFI- 
CATION FUNCTIONS  FOR  CATEGORIES  3  AND  4 
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SEASON  &  CATEGORY 

FOG  PERCENTAGE 

NO-FOG  PERCENTAGE 

DRY   (1+2) 
WET  (1+2) 
DRY   (3+4) 
WET  (3+4) 

85.1 
84.6 
71.0 
75.8 

82.1 
82.4 
75.5 
70.7 

WEIGHTED  AVERAGE 

74.9 

75.5 

TABLE  XXII.    QUILLAYUTE  FOG  AND  NO-FOG  DISCRIMINATION 
PERCENTAGES 
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CATEGORY 

NU2-IBER    &    CLASSIFICATION 

1 
2 

35    FOGGERS 
350    NO-FOGGERS 

3 
4 

20    FOGGERS 

97     NO-FOGGERS 

TABLE  XXIII.    SAN  DIEGO  DRY  SEASON  FOG  AND  NO-FOG  CASES 
FOR  THE  TEST  DATA  SAMPLE 
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CATEGORY 

NUMBER  &  CLASSIFICATION 

1 
2 

40  FOGGERS 
250  NO-FOGGERS 

3 

4 

76  FOGGERS 
150  NO-FOGGERS 

TABLE  XXIV.    SAN  DIEGO  WET  SEASON  FOG  AND  NO-FOG  CASES 
FOR  THE  TEST  DATA  SAMPLE 
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SEASON  &  CATEGORY 

FOG  PERCENTAGE 

NO-FOG  PERCENTAGE 

DRY   (1+2) 

51.4 

78.6 

WET   (1+2) 

70.0 

73.2 

DRY   (3+4) 

50.0 

72.  2 

WET   (3+4) 

80.3 

67.3 

WEIGHTED  AVERAGE 

68.4 

74  .3 

TABLE  XXV. 


SAN  DIEGO  FOG  AND  NO-FOG  DISCRIMINATION 
PERCENTAGES  FOR  THE  TEST  DATA  SAMPLE 
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APPENDIX  A 


A  GENERALIZED  FORTRAN  PROGRAM  FOR  APPLICATION 
OF  THE  CLASSIFICATION  FUNCTIONS 


C        VAR(I)  =  VARIABLE  RAW  VALUE 

C        N  =  NUMBER  OF  TERMS  IN  THE  CLASSIFICATION  EQUATION 

C        CF<I)  =  FOG  FUNCTION  COEFFICIENT 

C        CNF (I)  =  NO-FOG  FUNCTION  COEFFICIENT 

C        DF  =  FOG  DISCRIMINANT  SCORE 

C        DNF  =  NO-FOG  DISCRIMINANT  SCORE 

C        PPF  =  POSTERIOR  PROBABILITY  OF  FOG 

C        PPNF  =  POSTERIOR  PROBABILITY  OF  NO-FOG 

DIMENSION  VAR(N)  ,CN  (N)  ,CNF  (N) 

DF  =  CONST 

DNG  =  CONST 
C        READ  IN  THE  VARIABLE  RAW  VALUES 

READ(5,5)   (VAR(I) ,I=1,N) 
5  FORMAT (8F 10. 2) 
C        MULTIPLY  THE  N  VALUES  OF  THE  VARIABLES  BY  THEIR 
C        ASSOCIATED  COEFFICIENTS  AND  SUM 

DO  10  1=1, N 

DF  =  DF  +  CF(I)*VAR(I) 

DNF  =  DNF  +  CNF (I) *VAR (I ) 
10  CONTINUE 
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C        TO  AVOID  TRYING  TO  EXPONENTIATE  A  NUMBER  TOO  LARGE 
C        TO  BE  REPRESENTED  BY  THE  IBM  360  COMPUTER  DO  THE 
C        FOLLOWING: 

DMIN  =  AMINl (DF,DNF) 

DIFF  =  DF  -  DMIN 

IF (DIFF. GE. -174 .67)  GO  TO  20 

DIFF  =  0.0 
20  PF  =  EXP  (DIFF) 

DIFF  =  DHF  -  DMIN 

IF (DIFF. GE. -174 . 67)  GO  TO  30 

DIFF  =  0.0 
30  PNF  =  EXP  (DIFF) 

SP  =  PF  +  PNF 
C        NOW  FORM  THE  POSTERIOR  PROBABILITIES 

PPF  =  PF/SP 

PPNF  =  1.0  -  PPF 
C        NOW  OUTPUT  THE  PROBABILITIES  IN  ANY  DESIRED  FORMAT 

STOP 

END 
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